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TECHNICAL MEMORANDUM X- 64633 


A METHOD FOR NONLINEAR EXPONENTIAL 
REGRESSION ANALYSIS " 


INTRODUCTION 


The investigation of physical processes frequently requires the use of 
models that simulate or describe the processes. A model is often chosen so 
that certain variables interact in the model according to physical theories 
associated with the particular process. Formulation of a model often results 
in the form referred to as mathematical models. This is the familiar repre- 
sentation of the physical process by one or more equations that encompass the 
physical theory. A model equation contains identified independent variables 
and unknown parameters. Regression analysis is the statistical tool used to 
determine these unknown parameters, thereby providing an analytical repre- 
sentation of the experimental data. 

The general procedure in regression analysis is to take partial deriva- 
tives of a specific model-dependent minimizing function. These partial deriva- 
tives are taken with respect to each of the unknown model parameters. If the 
set of equations obtained by setting these partial derivatives equal to zero can 
be solved by the usual algebraic methods, the fitting or analytical representa- 
tion is accomplished. However, if these equations are transcendental in one 
or more of the unknown parameters , they cannot be solved by the usual alge- 
braic methods. 

The processes of particular interest in this report are those that can 
be described by decaying exponential forms. A mathematical model that con- 
tains more than one exponential term results in a set of transcendental normal 
equations if conventional forms of regression analysis are used. Thus, one 
usually resorts to iterative methods that require initial estimates for the 
parameters. The method described herein involves the least squares proce- 
dure, whereby the noplinear problem is linearized by expanding in a Taylor 
series. In this iterative method, we first develop a starting nominal guess 
for the model parameters. A correction matrix is derived and then applied to 
the nominal guess to produce an improved set of model parameters. This pro- 
cedure is continued until some predetermined criterion is satisfied. The 
number of iterations necessary for convergence is closely related to this 
criterion, the initial estimates, and the form of the exponential model. 



Additional information on the various methods of curve-fitting decay- 
type data to a sum of exponentials is given in References i through 5. Pro- 
cedures for obtaining the initial parameter estimates are discussed in 
References 6 through 8. It is noted that the initial estimate procedure herein 
is not restricted to equally spaced data. 

Application of the procedure is illustrated with data obtained from a 
particular process concerning the anodic oxidation of metals. In this process 
one expects an exponential or logarithmic behavior. From an analysis of the 
results, it is concluded that an adequate two-term exponential representation 
of the data is obtained. Thus , the analytical representation of the physical 
process data is accomplished using an exponential decay-type model. 


MATHEMATICAL THEORY FOR EXPONENTIAL 
REGRESSION ANALYSIS 


In general, we are given the set of observed values | (tj, fj 0 ) , 

(t 2 , f 2 °) , . . . , (t n> f n °) j. We assu me that the function to be fitted to these 
data is of the following form: 

— B t 

f C = A ie -Blt + A 2 e“ B2t + . . . +A m e m + K = f (A, B, K, t) , (1) 

Q 

where f represents the calculated value of the response and K, A^ , B^(i = 

1,2,..., m) are the 2m + 1 parameters to be estimated. The independent 
variable is time t . We first consider the simple case m = 1 and K = 0 
which results in a conventional least squares solution for the unknowns Aj 
and Bj . We have the following model: 



( 2 ) 


By taking the natural log of both sides , 


In f C = In A t - Bjt 


(3) 


2 


Let 


C = In A 1 


and 


(4) 


Y = In f° 


Then, 


Y = C - Bjt 


(5) 


The least squares solutions for C and Bj then become 


B, = - 


n n n 

m Y. t.Y. - V t. Y Y. 
i=i 11 i=i 1 i=l 1 


n 


n 


m 




( 6 ) 


and 


n n n n 

2 ‘i* 2 Y i - 2 2 *,*, 

i=l i=l i=l i=l 


n 


n 


(7) 


m - lA 


For m > 1 and K + 0 we proceed as follows. We write the parameters as 
initial approximations or nominal values plus unknown corrections; that is, 
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A t = A t + AAj 
B^ = B^ + ABj 

A 2 = A 2 + aa 2 
b 2 = b 2 + ab 2 


A 

= A 

+ AA 

m 

m 

m 

B 

r>~> 

= B 

+ AB 

m 

m 

m 


K - K + AK 



( 8 ) 
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For any assumed functional form, the following condition equation can be 
written: 



(9) 


where 

F = observed values of the response variable, 
f. = calculated values of the response variable, 

and 

V f . = residuals associated with the response variable. 

If we substitute equation (8) into equation (1) and the result into equation (9) , 
we obtain 


f°-f. (Aj + AAj, Bj + ABj A +AA , B + AB , 

1 1 \ m mm m 

K + AK, t) - V = 0 , (10) 


4 


or 


f. \ Aj + AAj , Bj + ABj , . . . , A + 
1 \ m 


AA , B +AB , 
mm m 


K 


+ AK, t ) = f° - V 
' l fi 


( 11 ) 


Expanding the left side of equation (11) in a Taylor series about the estimates 
Aj i Bj , , , , , A^ > B^ , I 
we have (i - 1, 2, . . . , n) 


Aj , Bj , . . . , A^ , B^ , K and neglecting higher order terms than the first, 


V_. = N. - F, . AA _ F„. AB 4 - F„. AA„ - F .. AB„ 

fi i li 1 2i 1 3i 2 4i 2 


F . . AA - F„ . AB - F„ , . AK 
2m- l,i m 2m, i m 2m+l,i 


( 12 ) 


where 


Sf. 


F„ = 


li 3A 


Of. 


F.. = 


4i dB, 


0f. 


F„. = 


2i 3B , 


3f. 


F_. = 


3i 3A 


3f. 


2m-l,i 3A 


m 


3f. 


2m , i 3B 


m 


(13) 


3f 

, = 1 

2m+ 1 , i 3K 
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and 


N. =f.° - f. (A, B, K, t) 
11 1 


( 14 ) 


The n equations given by equation (12) are the linearized condition or residual 
equations. According to the Gauss least squares principle, the best represen- 
tation of the data is that which makes the weighted sum of the squares of the 
residuals a minimum. Thus, the minimizing function is 

S =f (AA, , AB, , . . . , AA , AB , AK^j 
\ 1 1 m m / 

2 2 2 

= WV+WV+ +WV (15) 

1 fl 2 f2 ’ ' • n fn * U ' 


The 2m + 1 linear algebraic equations for determining the A increments to 
the initial estimates are now obtained by taking the partial derivative of S with 
respect to each of the unknown corrections and setting the result equal to zero. 
These equations, frequently referred to as the normal equations, are given by 


as 


= 2W ' V 


av av 

fl f2 

+ 2W V„ — — + 


3AA . 1 fl 9AA , 2 f2 3AA 


av. 

. + 2W — -f- = o 
n 3AA , 


as 


= 2W , V 


9V fl aV f2 

+ 2W V. — 1 — + 


3AB 1 1 fl 9AB 1 2 f2 3AB 


av. 

•• + 2 W n^t = ° 


(16) 


av.. av av. 

= 2W ,V., + 2WV-^ +. . . + 2W. ^ - 0 


3AK 1 fl 3AK 2 f2 3AK 


n 3AK 
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or 


7 



AA t 

ABj 

aa 2 


A 

[ (2m+l)x i] 



AA 

m 

AB 

m 

AK 


and 


Ni 


N 

[nxl] 


n 2 


N 

n 


( 21 ) 


( 22 ) 


By using the denotations for V and B as given by equations (19) and (20), 
respectively, we can also rewrite equation (17) as 

— T 

B W V = 0 , (23) 


where 




0 


0 


w = 

[nxn] 


0 



0 


0 


2 . 2 
ct /cr 
o fn 


( 24 ) 
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Solving equation (18) for V , 


V = N-B _ A 7 ” 


( 25 ) 


Substituting equation (25) into equation (23) , 


B T W(N-BA)=0 . (26) 


Solving for A , 

B T W N - B T W B A = 0 (27) 


or 


A = (B T W B)~ 4 (B T W N) 


(28) 


An improved set of values for the parameters is then given by equation ( 8) . 
The process or cycle is repeated to produce the corrections resulting from the 
second cycle. These corrections are then added to the estimates from the 
first cycle: 


A/ = AAj 1 + Aj 

B/ = ABj 1 + Bj 
A-2 = A A 2 + A 2 
= ^B2 1 + B 2 

A 1 = AA 1 + A 
m mm 

B 1 = AB 1 + B 
m mm 

K 1 = AK 1 + K 


(29) 
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These values represent an improved set of estimates to use for the third cycle. 
An iterative procedure is thus set up for improving the parameter estimates to 
any prescribed degree of accuracy consistent with the accuracy of the observed 
data. 


The standard algorithm is based on obtaining a nominal solution that, 
hopefully, converges to the correct solution. The algorithm is summarized 
as follows: 

~k 

1. Let A denote the kth nominal; linearize about A . 

2. Solve the resulting linear least squares problem. 

3. Use the new solution as the new nominal. 

4. Check for convergence. If convergence has not occurred, repeat 
steps 1 through 3. 

The standard deviation of each of the converged parameters is calculated from 



(30) 


(31) 


and 


1 cn 




[ (2m+l)xl] 


(32) 


41 


2m+l , 2m+l 
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The c elements in equation (32) refer to diagonal elements in the inverse 
matrix (B^WB)" 1 . 

PROCEDURE FOR OBTAINING THE INITIAL PARAMETER ESTIMATES 

"Peeling-Off" Approach 

An iterative method for nonlinear exponential regression analysis was 
developed in the previous section. Inherent in this method is a requirement 
for initial estimates of the parameters. This section presents a least squares 
"peeling-off" procedure for arriving at these initial estimates. 

Our assumed exponential model is of the form given by equation ( 1) . 

Generally speaking, if we plot decay-type data in the form In f° against t 

where f° is the observed response and t is the independent variable, then 
for large t the curve is approximately a straight line. Consider the following 
Figure 1. 



Figure 1. Logarithmic time decay illustration. 

If we fit a straight line to the last three data points by the method of least 
squares, the assumed form is 


In f = - B t + D 
m 


However, this is equivalent to the equation 


-B t 

f° = A e m 
m 


(33) 
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where 


D = In A 


m 


(34) 


The coefficients B and D are given by 
m 


n-2 n-2 n-2 n-2 

1 1 1 Z “n 0 - Z t, z 


D 


i=n i=n 


i=n i=n 


n-2 /n-2 

3 z *i - z 


i=n 


i=n 


(35) 


and 


B = - 

m 


n-2 


n-2 n-2 


3 Z \ ln f i° - Z fc i Z ln f i° 


i=n 


n-2 


l-n i=n 
'n-2 


3 Z h 2 - Z *. 


i=n 


i=n 


(36) 


Thus, we have determined the least squares values for A and B . We 
’ mm 

can then obtain values for the residuals from (i = 0,1,2) 


-B t . 

„ r o ~ m n-i 

R . = f . - A e 
n-i n-i m 


(37) 


We now take the next three data points at t . , t . , and t _ and subtract 

n-3 n-4 n-5 

-B t 

~ m ~ o 

the corresponding term A^e and the arbitrary constant K from f to 

obtain the following residuals 
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o o 


O O 


R = f° - A e 
n-3 n-3 m 


-B t „ 
m n-3 


- K 


-B t 

R = f» - A e mn - 4 - K 
n-4 n-4 m 


-B t _ 

R . f» - A e m I1 - 5 - K 
n-5 n-5 m 


( 38 ) 


Next, a straight line is fitted to the data for In |R^ | against t^ . , i - 

3,4,5. This determines A , and B , . The residuals for the next 

m-1 m-1 

three data points are determined from 


R = f 


- fO 


n-6 n-6 m-1 


-B .t -B t 

~ m-i n-6 ~ m n-6 ~ 

A ,e - A e - K 


m 


-B t _ -B t _ 

R = f 0 - A e m - ln - 7 -A e m D_7 - K 

n-7 n-7 m-1 m 


(39) 


R 0 = f° . - A ,e 
n-8 n-8 m-1 


-B .t -B t „ 

m-1 n-8 ~ m n-8 ~ 

- A e - K 

m 


We now fit a straight line to In J against t^ . , i = 6,7, 8 . This 


determines A _ and B „ . We continue this process until all the A's 
m-2 m-2 


_ />*/ rs . » 

and B's are determined. In general, Aj and Bj are determined from a 
set that contains more than three data points. That is, the points remaining 

r*j r*J 

after A 2 and B 2 have been determined are used for determining A t and Bj . 
It is noted that since three points are chosen as a minimum for calculating a 

particular A and B , we must have n S: 3m where n is the number of data 
points. At this point we calculate the weighted sum of squares of the devia- 
tions using the initial parameter estimates 
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(40) 


Fw = 


| 1 w ‘( f *°- , *°) 


where 


W. = i/crj 

1 fi 


f = observed response 


and 


f = fitted response. 


Iteration Philosophy 

The iteration logic can be summarized by the following steps: 

1. Repeat the process, but use the last four data points to obtain the 

initial estimates of A and B . The next three data points are used to 

mm 

obtain A , B . , etc. , for other A’s and B's . This yields a second 
m-i m-1 ~ 

set of parameter estimates from which we can calculate another F ; call it 

f 2 . 

2. Use the last five data points and obtain a third set of parameters 

r>j 

which yield F 3 . 

3. We continue this, always keeping three points as a minimum in 
determining A^ and B. . 

4. Increase the constant to K + 0. 05 K and repeat steps 1 through 3. 
Terminate when the constant reaches a predetermined value. 

5. Parameter estimates yielding the minimum F are chosen as the 
initial estimates. 
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COMPUTER PROGRAM DEVELOPMENT 


The com puter programs to implem en t the previously developed theory _ 

for exponential regression analysis were organized and developed according to 
two general types of exponential models. One concerns a single exponential 
and the sum of exponentials without a constant, and the other concerns the sum 
of exponentials with a constant included. Two highly flexible computer pro- 
grams were thus developed for the MSFC UNIVAC 1108 digital computer. Each 
program contains double-precision capability and SC-4020 plotting procedures. 
In addition, the "peeling-off" procedure for obtaining initial parameter esti- 
mates is an integral part of each program segment. 

The logic flow for Programs I and II is depicted in Figures 2 and 3. 

The parameter NCASES is the number of cases of data processed in each pro- 
gram. The models that can be investigated in Program I are 


Model I: f C = Aie -2 ^ + A 2 e _B2t + Age -03 * + K 

and 

Model II: f C = A ie ~ Blt + A 2 e~ B2t + K 


Those models that can be investigated in Program II are 


Model III: f C = A ie ' Blt + A 2 e Bzt + A 3 e~ B3t 


Model IV: f C = A ie Bjt + A?e~ Bzt 


and 


Model V: f C = A ie " Blt 


The characteristics of both programs are summarized in Table 1. It should 
be noted that both programs can be easily extended to include additional expo- 
nential terms if desired. 
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Initial Estimates 
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Figure 2. Program I block diagram summary. 
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Figure 3. Program II block diagram summary. 






















TABLE 1. COMPUTER PROGRAM SUMMARY FOR 
EXPONENTIAL CURVE FITTING 


Program 

Number & 

Exponential Model Fitted 

I 

.(1) A. , B. , K : i = 1,2,3 

l l 


(2) A. , B. , K : i = 1,2 

i l 

II 

(1) A. , B. : i = 1,2,3 


(2) A. , B. : i = 1,2 


(3) Aj , Bj 


a. Each program has a built-in capability for 
obtaining initial parameter estimates. 


The parameter ITERM is used to determine control transfer in each 
program. For example, if ITERM has the value one, parameters for the 
particular model associated with it are determined. A value of zero indicates 
that parameters for the associated model are not determined. A double- 
precision matrix inverse routine using the Gaussian elimination procedure is 
used in each program. 


RESULTS AND CONCLUSIONS 
Discussion of the Physical Process 

In an open circuit transient analysis of the anodic oxidation of metals, 
one expects an exponential or logarithmic behavior. This fact is evident from 
the experimental data when the layer is passive; that is, normal growth is 
taking place. In some cases, however, growth is truncated by the onset of 
-oxygen evolution. In this region we have electronic conduction in addition to 
a small amount of anodic or growth conduction. During the open circuit break, 
therefore , we expect the conduction to be initially dominated by electrons and 
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as the voltage decreases , the conduction becomes primarily anodic conduction. 
Thus , the transient analysis has to have the capability to consider two or more 
conduction mechanisms, each with different relaxation times. 

Basically, the voltage data from the anodic oxidation process are classi- 
fied into four sets. Each set of data represents a different time held in oxygen 
evolution. When plotted as a function of time, these voltage data exhibit an 
exponential 'or decay-type behavior. The response data are represented as an 

o 

observed voltage, which is designated as f . 


Discussion of Results 

The observed decay data were processed through both Programs I and 
II to assess the validity of the various exponential models. A summary of the 
initial estimates for the various assumed models is presented in Tables 2, 3, 

and 4. The minimum value for the F quantity associated with the selected 
initial estimates is given in these tables. 


TABLE 2. INITIAL ESTIMATE SUMMARY FOR TWO-TERM 
EXPONENTIAL MODEL PLUS A CONSTANT 



Set 1 

Set 2 

Set 3 

Set 4 

Groups of Initial 
Estimates 

441 

441 

441 

441 

Group Number with 

49 

44 

61 

63 

Smallest F 





Smallest F 
Value 

0. 2083 

0. 3113 

2. 9021 

5. 986 
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TABLE 3. INITIAL ESTIMATE SUMMARY FOR TWO-TERM 

EXPONENTIAL MODEL 



Set 1 

Set 2 

Set 3 

Set 4 

Groups of Initial 
Estimates 

63 

63 

63 

63 

Group Number with 

52 

63 

61 

59 

Smallest F 





Smallest F 
Value 

0.0827 

0. 0455 

0. 01185 

0. 02952 

TABLE 4. INITIAL ESTIMATE SUMMARY FOR THREE- 
EXPONENTIAL MODEL PLUS A CONSTANT 

TERM 


Set 1 

Set 2 

Set 3 

Set 4 

Groups of Initial 
Estimates 

420 

420 

420 

— 

Group Number with 

59 

47 

14 

— 

Smallest F 





Smallest F 
Value 

0. 2132 

0. 3532 

57.45 



Numerical problems were encountered for both the three-term model 
and the three-term plus a constant model. As indicated in Table 5, divergence 
occurred for the set 2 and set 3 data. It is noted that the determinant of the 
coefficients for solving for the correction matrix was 0. 911 x 10 -18 in one 
case and 0. 162 x 10 -2 in the other case. Convergence failed to occur for the 
set 2 data even with relatively good initial estimates. 
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TABLE 5. RESULTS FOR THREE-TERM EXPONENTIAL MODEL PLUS A CONSTANT 


Xi 

T3 <U 
<D +} 

O g 

^ -5 

Q, •+■» 

















Set - 

« w 









Initial 

Estimate 
















. 

Set 3 

Improved 

Estimate 21 








CO 

TH 

1 

o 

tH 








tH 

tH 

05 

o 

Initial 

Estimate 

0. 289545 

0. 006907 

0. 237985 

0. 000905 

0. 777985 

086900 *0 

0. 21350 

0. 9705 

Set 2 

Improved 

Estimate 3 - 








i 50 

1 ^ 








§ . 

o' 

Initial 

Estimate 

0. 524024 

0. 121633 

0. 036133 

0. 065537 

0. 473872 

0. 011831 

0. 0770 

0. 0761 

Set 1 

Improved 

Estimate 

0. 550958 
0. 0319 

1. 003510 
0. 3064 

0. 571808 
0. 0200 

0. 063911 
0. 0070 

0. 108318 
0. 0066 

5000 *0 
U68I0 *0 

05 05 

00 co 

CO O 

o o 
o o 

0. 0321 
7. 90 

Initial 

Estimate 

0. 506027 

1. 144510 

0. 525248 

0. 117757 

0. 231458 

00 

co 

° 1 
1 1 
o 

0. 0525 

0. 0591 

Data Set 
Number 

Parameter 

and 

Error Estimate 

rH 

< b 

ffl 1 b 

C\J 

< b 

CVJ 

M b 

CO 

<: b 

CO 

pq b 

K 

a 

°f 

Determinant 

Value 


21 


Diverged. 

Numerical problems for assumed model. 



The data in Table 6 summarize the results for the two-term plus a 
constant model. As shown in this table convergence failed to occur for the 
set 3 and set 4 data. The initial estimates for the set 1 and set 2 data showed 
some disagreement with the cycle 1 estimates. The initial estimate values 

for set 1 and set 2 are significantly smaller than the cr^ values at the end of 

cycle 1. Rather high error estimates for the A t and parameters are 
also evident in this table. 

Results obtained for the two-term exponential model are perhaps the 
most encouraging from the standpoint of adequately describing the data. As 
indicated in Table 7, highly accurate parameter estimates were obtained. The 
converged estimates represent an improvement over the initial estimates with 
the exception of the set 4 data. Here it is noted that = 0. 0215 for the 

initial estimates as compared to = 0. 0600 for the improved estimates. The 

initial estimates were thus chosen as the representation for the set 4 data. 

The models that appear to adequately describe the observation data are 

Set 1: 

Set 2: 

Set 3: 

and 

Set 4: 


f c = 0. 924798e-°- 1607831 + 0. MOM*-*- 0107511 


f C = 0.629470e-°' 230231 + 0.676240e-°- 016421 


f C = 0. 598702e°‘ 1263261 + 0. 707297e-°' 0061321 


f C = 0. 424970e-°- 1093681 + 0. 897037e-°- 0058761 


The observation data and the models evaluated at the corresponding time point 
are presented in graphical form in Figures 4 through 7. Residual data are 
presented in Figures 8 through 11. These figures indicate an adequate model 
representation of the data. It is concluded that changing the exponential func- 
tional form for these data to one other than a two-term model is not warranted 
in view of the problems encountered with the other models. 

The coefficients for each model are plotted in Figure 12 as a function 
of the time to oxygen evolution associated with each set. These data enable 
one to simulate the physical process using a two-term exponential model at 
conditions other than those tested. 
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TABLE 6. RESULTS FOR TWO-TERM EXPONENTIAL MODEL PLUS A CONSTANT 
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TABLE 7. RESULTS FOR TWO-TERM EXPONENTIAL MODEL 
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Figure 4. Observed and computed response for Model IV analysis 

on set 1 data. 
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TIME 


Figure 5. Observed and computed response for Model IV analysis 

on set 2 data. 
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Figure 6. Observed and computed response for Model IV analysis 

on set 3 data. 
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RESPONSE 



TIME 


Figure 7. Observed and computed response for Model IV analysis 

on set 4 data. 
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TIME 


Figure 8. Residuals for Model IV analysis on set 1 data. 
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Figure 9. Residuals for Model IV analysis on set 2 data. 
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Figure 10. Residuals for Model IV analysis on set 3 data. 
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TIME 

Figure ii. Residuals for Model IV analysis on set 4 data. 
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APPENDIX 

COMPUTER PROGRAM DOCUMENTATION 


This appendix presents operational information on the UNIVAC 1108 
computer programs concerning the regression analysis application of expo- 
nential models to decay-type data. The general organization of the operational 
version of the two developed programs is depicted in Figure A-l. Since both 
programs are basically similar, only information concerning Program II is 
presented. A complete program listing, job card example (Fig. A-2) , input 
preparation, and sample output are included. 



Figure A-l. Program organization. 
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1108 RUN REQUEST & INSTRUCTIONS 


NAME (LAST & INITIAL) 

JUNKIN. B. 

OP 

w 

JOB« 

320600 

PROD # 

BJN* 

313 

BLOCH 

4663 

RM# 

A-217 

RUN-ID 

EXP REG 

RUN 

1 OF 1 

EST. CPU RUN TIME 
-Q_ HRS. 3 _ MINS. 

CORE SIZE 

24K 

PUNCH 

N/A 

i 

C 

3 COMPILE 
3 EXECUTE 
] SORT 

□ 

(B EXEC VIII 

LANGUAGE 

FOR4 

MAX. 

PAGES 

1,500 

□ SPECIAL FORMS 
TYPE COPIES 


DOES THIS JOB HAVE A RESTART PROCEDURE? IDYES D NO 



Set 1 Data: 2 exponentials 

□ OVER 


MICttO FILM 

COPIES 

COPY 

FLO 

OPER. IN IT. 

IFILES 

^FRAMES 

P 

F 

P 

F 

SBQ.lt 









OPERATOR COMMENTS! □ SEE TECH. □ SEE OPER. 


MSFC - Form 3019 (**v Augu.t 1909) 


D OVER 


Figure A-2. Job card example 



Description of Data Deck Input Parameters 

The card immediately following the @XQT card is the first card of the 
input data deck. This card specifies the number of cases (NCASES) of data 
that are to be processed. The format is of the form +XX(13) and appears in 
columns 1 through 3. The information between the $ INPUT card and the 
$ card is associated with a specific set of data and is input under the non- 
executable NAMELIST statement. For example, the input statement in the 
program is 

NAME LIST/INPUT/T , Y , NN , TL , TR , YB , YT , YB 1 , YT i , VARY , TOLER , 

ITERM4, ITERM3, ITERM2. 

The forms that the input data take include variable name and sub- 
scripted variable. In the usage above , T and Y are subscripted arrays and 
the remaining variables are simple variable names. The specific format of 
the data can be either integer constants (i. e. , +218) or real constants (i. e. , 

1. 85921E+00, with or without the E notation). The description of the vari- 
ables in the NAMELIST statement follows. 


T 

Y 

NN 

TL 

TR 

YB 

YT 

YB1 

YT1 

VARY 

TOLER 


— array containing the values for the independent variable 

— array containing the values for the dependent variable 

— number of data points 

— left plot limit for the horizontal T axis 

— right plot limit for the horizontal T axis 

— bottom plot limit for the vertical Y axis 

— top plot limit for the vertical Y axis 

— bottom plot limit for vertical residual axis 

— top plot limit for vertical residual axis 

— cr £ , variance for dependent variable 

— iteration parameter 
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ITERM4 

ITERM3 

ITERM2 


— control parameter for Model III 

— control parameter for Model IV 

— control parameter for Model V 



Computer Listing of Program II 
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00 

COO- 

p 

L 1 3 1 1 6 

DEL 

Lean 

0 

C23&I4 

OLTEr 

-oco 

ft 

C-2 3C7 6 

OF 

COCO 

0 

L2 1 642 

01 DEN 

0300 

R 

-23 1 33 

OTOL 

C Go. 

b 

-IC.Gl 

01 

204)6 

0 

0*2574 

02 

-003 

U 

012612 

03 

ooco 

0 

j2 1 5G- 

D 3 3 

OOOO 

ft 

-1 3754 

F 

GOG- 

1 

y23b4* 

1 

vzub 

1 

C23-4S 

10 

cGGO 


023134 

INPUT 

0LC0 

I 

i.2304 1 

1 TEftM2 

ouoo 

I 

J 2 3 C 4 0 

I I ERM 3 

OL Ui/ 

I 

£23.37 

1 TLRMft 

0606 

1 

023l 17 

I " 

-OoO 

1 

02311 1 

i X 

GOOD 

1 

-23 1 1 2 

I 2 

□uuo 

I 

£23-74 

J 

bti OiS 

1 

C23IB* 

JB 

Co-6'0 

I 

023 t 5o 

JO 

CC-LG 

I 

023C73 

JJ 1 

00 CO 

I 

C 2 3 ijS 6 

JN 

OUOO 

I 

323 10Q 

JSEL2 

GCU- 

1 

.23.77 

JSEL3 

-060 

i 

.23.46 

Jl 

00.3 

1 

C 2 3q55 

J2 

COCO 

1 

L 2 3 — 75 

A 

00-0 

I 

-23101 

ACYCLt 

GCuj 

1 

u231(]0 

Vk 

1006 

1 

023.7 1 

KICK 

.000 

1 

02311 3 

K 1 

OOCC 

I 

0231 1 4 

A 2 

0003 

1 

0231 1 5 

K3 

ljoc 

I 

023116 

K 4 

00 —O' 

1 

0231 1 c 

n 

OOwO 

1 

023 1 q7 

M 

GOOC 

I 

-23o44 

NCASES 

OUOO 

1 

-23-72 

NKOT 

COO- 

1 

- 2 3 l 26 

Nil 

souo 

1 

023.42 

p 

£OJ0 

u 

013134 

ftES 

00 CO 

ft 

-22160 

RSS 

ouco 

H 

J231 24 

5b A 1 

GO 0- 

K 

■123126 

SGA? 

L u Uu 

K 

.23131 

56*3 

COJO 

ft 

023125 

5GB 1 

ouo 

ft 

tZ3 1 27 

SGB2 

OUOO 

ft 

J23 1 3 1 

S b B 3 

(‘Cub 

R 

0231 C‘4 

S lb Y IN 

wC60 

K 

023121 

sun 

£0-0 

ft 

023123 

S1GY2 

ccoc 

G 

-23016 

5H 1 

cuoo 

D 

023-20 

SM2 

oco- 

D 

.23.22 

SM3 

60-6 

D 

023.24 

SH4 

jO.O 

ft 

G2307G 

5UM0EV 

0000 

ft 

-23103 

SUMS- 

OJCO 

ft 

323120 

SUM 1 

Cbtov 

K 

123122 

SUM2 

wOClO 

0 

GOOoCC 

T 

30C0 

ft 

021334 

TEHP 

GliiO 

ft 

-23-27 

XL 

oooo 

R 

£23-36 

TOLEH 

IGU- 

ft 

123132 

10L1 

.CuC 

ft 

023.3,. 

1ft 

<3C jO 

ft 

022324- 

TS 

CL CO 

ft 

-23-35 

VARY 

OuuO 

R 

-0C620 

A 

COO- 

0 

tvt 76** 

X 1 

CL-- 

0 

0*01274 

X 2 

£030 

U 

CC 1 6 0 4 

X3 

OLOC 

D 

-CQ31G 

Y 

oooo 

ft 

-23C31 

YB 

OLC- 

N 

v^3j3 J 

rb i 

OcCw 

0 

013444 

tc 

00. C 

ft 

022634. 

res 

3000 

ft 

-2247C 

YS 

OJCO 

R 

-2 3-32 

Y T 


CCOu R Z-2 3-34 r T I 


LG ICC 

i • 

C 


-OlOfc 

2* 

c 

NUMERICAL TEChMUOE FOR EXPONENTmt REGRESSION ANALYSIS 

CCI-t 

3 • 

c 


toll 1 C 1 

4* 


01 MENS ION T* l.C* • Y f 1 9; 1 ,*«|0C> .X 1 1 1CU) .X2| ICO) i A3I ICO) >B( lOj .7 1 

.51-1 

5* 


-Bll 7 , |£v) ,8NI 1 -G, 1 1 .OH 100,7 1 ,021 7 , 1 1 ,D3I 7 .7 ) ,D 1 4 9 I , DELI 7 , | ) , 

i.ClC-1 

6 • 


-RES< !;£•> »YC 1 I.CI ,F 1 35C» ,AR« 7 , 350) , TEhP < ( ,';C ) . 0 33 1 7 , 7 > , 0 1 OE N I 7 , 7 ) 

Gild 

7 • 


-BCO T 1 1 2 1 ,BCO V 1 1 2) .BCOKI 1 2) ,DC( 36 ) 

0CIC-3 

8* 


DIMENSION RSSl l£j) ,1S1 IOC) , YS1 ICC) ,YCS1 ICC) 

-01-4 

V. 


DOUBLE PfttCiSiON X | , x 2 , X 3 , Y f . «E 5 , A l | , 8 | |, A Z2 , 0 22 , A 3 3 , 83 3 » DE L 

-cits 

lw* 


DOUBLE PRECISION B f B T , BN , 0 1 • D2 , 0 3 ,0 »0D , DE T ER . 0 3 3 , D 1 OEN , T , Y 

-.1-6 

ll« 


DOUBLE PRECISION 5M 1 , SM2 ,SM3 ,SM4 

CC 1 1)7 

!*• 


NA*ELIST/1NPU!/T»Y.Nn,TL,1R.Yb.YT,YBI ,YT1 , VARY .TOLER. ITERM4, 

iiU7 

13-' 


*• 1 T £ RM 3 » 1 TLKM2 

tone 

14- 


CALL 1DENT1935) 

com 

lS- 


INTEGER P 


44 



Sit) 12 

16* 

1 “V 

FORMAT 

I 1H1 1 

Ou 113 

17. 

1 Vi 

FORMAT 

( // ) 

tans 

1 6. 

14.2 

format 

(83 H N0N-LIn£A 8 EXPONENTIAL REGRESSION ANALTSlS USING AN IT 

ijlH 

1 9. 


-E. ft A T 1 VE 

CORRECTION PROCEDURE) 

ub 1 1 5 

20. 

U-3 

format 

( 36H ASSUMED MODEL [S THREE EXPONENTIALS) 

aa i it 

2)« 

lvH 

format 

( 33H INITIAL ESTIMATES FOR PARAMETERS) 

Cut 1 7 

22* 

I c 

format 

< 3X , 3mA 1 *E | J . 6,2x , J hB)=E 1 3, 6.2X, 3HA2-E 1 3. 6, 2X , 3H82*EI3Y6 , 2x 

cen7 — 

— 2 3« — 


-,'3HA3»E 13.G ,2X,3HB3.£13»6> 

60126 

21. 

106 

FORMAT 

( I5H CTCLE NUMBER » I 3 ) 

60121 

2b. 

1 - 7 

FORMAT 

( 3X i 22h PARAMETER CORRECTIONS) 

06 122 

26. 

168 

FORMAT 

( 3X ,1hOA I »E| 3.6 , | X. 1MDB 1 =E ( 3.6 , I X . 1hOA2=E| 3.6, |X , 1H0B2.E 1 3. 

66122 

27. 


”6. 1 X , 1 H D A 3 = E ) 3.6 , IX • 1 M OB 3* E 1 3.6) 

00123 

28. 

109 

format 

(3x,29h improved parameter estimatesi 

03121 

29. 

Do 

FORMAT 

( 3 X , 3 7m RESPONSE VARIABLE STANDARD DEVIATION) 

0CI25 

30. 

1 1 1 

F 0-R M A T 

( 3X , 6HS I G Y |*E | 3.6 .3X i6HSIGY2“E 13.6 ) 

66126 

31. 

112 

FORMAT 

I 3X , 3 1 H COEFFICIENT STANDARD DEVIATION) 

061 27 

32. 

1 13 

FORMAT 

t3X,lHSAI=El3,6,|X,lHSB|.E|3.6,|X,1HSA2«E|3.6, l X , 1 HSB 2 «E t 3 . 

C6127 

33* 


-6,|X,1HSA3=£|3.6,|X,1hSB3«E13.6) 

CCI3C 

31* 

1 11 

format 

( 3X , 6HT0LER.E 13.61 

66131 

35. 

1 lb 

FORMAT 

( 3 X , 5hDT0L*E 13.6) 

60132 

36. 

116 

format 

( 3X , 29HCONVERGENCE HAS BEEN ACHIEVED) 

CSI33 

37* 

117 

FORMAT 

( 3x i i Bh values in last cycle are final parameter values) 

CO) 39 

38. 

1 it) 

FORMA T 

( 131 

66) 35 

39. 

1 ) V 

forma t 

(siH initial estimates for three-term exponential mooed 

00 36 

1C* 

12c 

format 

<SX,lhI.8XilMF.11X,2HAl.l3X.2HBI,13X.2HA2,|3X,2HB2.l3x,2HA3 

0036 

1»« 


)3X,2HB3I 

6CI37 

12» 

1 2 | 

format 

(3X.I3.8(E|3* 6,2X1) 

66 1 ‘tO 

13. 

) 22 

FORMAT 

( 3x , 3 1 rt the initial estimates used are) 

60 *• 1 

11. 

123. 

FORMAT 

I 3X , 2m 1 * ) 

00112 

15. 

121 

FORMAT 

(19h initial estimates for TWO-TERM EXPONENTIAL MODEL) 

00l13 

16. 

125 

FORMAT 

lbX,|hl,BX,iHF,|1X,2MA|,l3X,2MB|,|3X,2HA2,13X,2HB2) 

00111 

17* 

1 2 6 

format 

( 3X, I 3.6(EI3-6,2X ) I 

66 1 15 

18. 

127 

format 

(3DH INVERSE TIMES ORIGINAL MATRIX) 

■60116 

19. 

128 

format 

(3X,7|E|3.6,IX)) 

66)17 

St* 

129 

format 

(3X.I8H DETERMINANT VALUE) 

6050 

SI. 

1 J'v 

FORMAT 

(31H ASSUMED MODEL | 5 TWO EXPONENTIALS) 

Ob) 5 1 

52* 

1 3 1 

FORMAT 

<3*i3HA|*E|3.6 t 3X,3HB|»E)3.6»3X,3HA2*E|3.6,3X,3HB2*El3.6) 

6052 

S3. 

1 32 

FORMAT 

<3A,6h0ELAI*£l3.6,3X,6HDELBI*E)3.6,3X,6H0tLA2'EI3.6,3X,6H0E 

6052 

SI* 


- L B 2 °E 1 3 

.6) 

06 1 5 3 

55* 

1 3 3 

format 

(3X,3H»|.E|3.6,3X,3H B l“E|3.6»3X 1 3HA2*E|3,6,3X,3HB2*E|3.6l 

66)51 

56. 

131 

FORMAT 

l3X,5HSGAl«E»3.6,3X,5HSGBI»£l3.6,3X,5HSGA2«ei3.6,3X,SHSGB2» 

6051 

S 7 • 


•E 1 3 • 6 ) 


00)55 

se. 

1 35 

format 

(3X.I0H RESIDUALS) 

60)56 

59. 

136 

FORMAT 

(3X,t|3.6) 

66) 57 

6u* 

137 

format 

( 1 8M I NVERSE uf bt.wb ) 

60 ) 6u 

6 1 * 

1 30 

format 

( 9 h B MATRIX) 

too 

62* 

139 

format 

(SIM OBSERVED response, COMPUTED RESPONSE and RESIDUALS) 

66)62 

63* 

Hw 

format 

16X,8m0BSEKV£D,6X,8H00MP0TED>6X,8HRESIDUAL) 

06) 63 

61. 

111 

FORMAT 

I2IH end OF CrcEE NUMBER I3> 

COli** 

65. 

112 

FORMAT 

(76H OBSERVED RESPONSE, COMPUTED RESPONSE, AND RESIDUALS US 

60)61 

66. 


-ING INITIAL ESTIMATES) 

60 1 6 5 

67. 

113 

FORMAT 

( 3X , 7hS I GT 1 N*E | 3.6) 

00)66 

68. 

111 

FORMAT 

1 7 H BT.V5.bl 

60)67 

69. 

lib 

format 

(33H ASSUMED MODEL IS ONE EXPONENTIAL) 

00170 

70* 

1 16 

format 

(2UH PARAMETER ESTIMATES) 

Cb)7> 

71* 

117 

format 

(3X,3ha|>EI3.6,2X,3hD)>EI3.6) 
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6617 2 

72* 


OAT A- 4 SCOT I 14 .1 *1 . 12) / 6HT1 ME , 11 *6H 

/ 

66174 

73* 


DATA (BCOTU) ,l*l*l2)/6HRESf'0N,6HSE 

, 1 C • 6 H 

Si Ol 7 6 

74* 


DATA <BC0K<I),I«l,l2)/6HR£SlDU,6HAL 

, 1 0»6h 


75* 


READ (5.118) NCASES 


C&203 

76* 

2oc 

T-RITE (6, ICC) 


60 205 

77* 


REAO (S. INPUT) 


002 1 U 

76* 


kRITe (6, IN Pul ) 


30213 

79* 


IF « JTERM2.Ne.64 GO To 706 


002 15 

8u* 


IF < lTERM4.Eo,o) GO TO 4000 


u u 2 1 5 

6 i * 

c 



00215 

8 2* 

C INITIAL ESTIMATES f OK THRft EXPONENTIALS 


00215 

83* 

c 

STRAIGHT LINE Fll TO LOG Y VS. TIME 


00215 

84* 

c 



00217 

85* 


I J = s; 


00220 

86* 


J 1 *2 


00221 

67* 


B J 1 =3. 


03222 

8 6* 

3 u CO 

I J= I J+ l 


63223 

89 • 


SMI*C. 


63224 

VC* 


5M2=C. 


002 25 

9 1 • 


SM3>*0. 


0022 6 

92* 


SM4=0. 


00227 

93* 


JJsNH-J 1 


00230 

94* 


DO 30C I = J J . N N 


00233 

95* 


SM1*T(1)«*2+SNJ 


60234 

96* 


SM2=DL0G ( Y ( I) ) +SM2 


00235 

97* 


SM3“T « I > +SM3 


00236 

98* 

3wc 

SM4-T < I > *DLOg< Y < 1)) ♦SM4 


CC240 

9 V* 


CAA3*BJI*SM)-SM3*»2 


002 4 1 

I0u» 


CA3*(SM|*SM2-SM3»5M4)/CAA3 


6024 1 

151 * 

c 



Q024 1 

102* 

C ESTIMATES FOR A3, 83 


00241 

1 C3* 

C 



W02H2 

1 64* 


B3*-((BJI*SMM-SM3*SM2)/CAA3) 


002^3 

105* 


A3=EXP ( C A3 ) 


00244 

1 i>6* 


J2«J| *3 


00244 

157* 

C 



0024 4 

108* 

C STRAIGHT LINE FIT TO LOG RESIDUAL VS. TIME 


00244 

l C-9 • 

C 



0024 5 

1 1 u* 


■J J*NN-J2 


002 4 6 

111* 


jN*NN- j | - 1 


00247 

l 12* 


SM1-C. 


06250 

1 1 3* 


SM2 = 0 . 


00251 

1 1 4* 


SM3*o. 


00252 

1 15* 


SM4=o. 


0025 3 

1 1 6* 


DO 3d 1= JO.JN 


UU256 

1 1 7 • 


X 3 < I ) «-83*T ( 1 , 


60257 

1 1 e* 


RES(I)»Y(I)-A3*D£XPlX3fI)) 


66260 

1 19* 


RES< I )*0A6S<KES< i ) I 


0326 1 

1 2c* 


SM I *T~< 14 **2*Sm | 


0u262 

12 1* 


SM2=0L0G (RES ( 1 ) ) +SM2 


63263 

122* 


SM3«T ( I ) *SM3 


50264 

123* 

3o! 

SM4*T(()*0L0G(RESlr , )+SM4 


63266 

124* 


C A A 2 a 3 • *5M 1 -SM 3 * *2 


60267 

125* 


Ca2*(SMI»SM2-SM3*SM4)/CAA2 


60267 

126* 

C 



66267 

127* 

C ESTIMATES FOR A 2 , 8 2 
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00267 
002 70 
0027 1 
0027 t 
002 7 1 

126* 
12 9* 
1 3c» 
l 3 i * 
132* 

C . 

B2*-t(3.*SH-4-SM3*SM2)/CAA2) 

A2«EXP<CA2> 

C 

C STRAIGHT LINE KIT TO LOG RES I DUAL VS. TIME — 

0027 1 
002 7 2 

13 3* 

1 3-4* 

-i 

JN«NN- J2-1 

002 7 3 

135* 


8 J 2 * J N 

002 7 H 

136* 


SMI = cj. 

00275 

137- 


SM2*c. 

00276 

136* 


5H3 = j. 

00277 

139* 


SMH = ij . 

3030C 

l-tg* 


DO 302 1= 1 , JN 

003 0 3 

» -4 » • 


X 3 1 1 ) = -B 3 * T t I ) 

0230-4 

1-42* 


X 2 4 I ) =-B2*T ( | ) 

00305 

1-43* 


RESm=Y(D-A2*DtXP<X2<n>-A3*DEXP(X3m) 

00336 

1-4-4* 


RES ( I >»DABS4RES( I ) ) 

O03O7 

1-45* 


SM1*T ( I » • » 2 + S M 1 

003 tO 

l *46* 


SM2=0L0Gl RES (1)1 *Sm2 

0031 i 

1 *47* 


S M 3 * T < I ) + S M 3 

003 1 2 

1 -48* 

3c2 

SM-4*T ( I > *DlUG 4 REStl » ) *SMR 

003 

1-49* 


CAAI=BJ2*SM1-SM3.*2 

C03 1 5 

ISu* 


CA l* ( SM 1 *541 2-5113 *SM*4 ) /C A A 1 

%»ti3 t 5 

15 i* 

C 


« j3 1 5 

152* 

C ESTIMATES FOR A i ,B I 

003 1 5 

153* 

C 


003 t 6 

15-4* 


B1=*-((BJ2*SMr-sm3*SM2)/CAA1> 

Ow3l 7 

155* 


A 1 = EXP ( CA 1 } 

CC317 

156* 

c 


003 1 7 

157* 

C WEIGHTED sum of squares OF DEVI AT IONS ,W = 1 . 

00317 

150* 

c 


00320 

159* 


DO 3^3 1 = 1 , N u 

00323 

1 6c* 


X 1 ( 1 ) =-B 1 *T ( I i 

0032*4 

16 1* 


X 2 ( 1 ) =-8 2 • T ( 1 1 

30325 

162* 

3u3 

X 3 ( I ) = — B 3 • T ( I ) 

00327 

163* 


SUMOEV=J. 

SU33U 

16-4* 


00 3cv 1=1 , Nn 

00333 

165* 


Bilill =OEXP 4 X i 4 I )> 

003 3-4 

166* 


B 4 I , 3 ) = D £ X P 4 X2 4 I )) 

00335 

167* 


B 4 I ,S ) aOtXP 4 X 3 4 I ) ) 

00336 

1 68* 

3o-4 

SUMDEV=(Y(I)-AI*B(I»1)-A2*B<J.3)-A3*8CI,5))»*2*SUHDEV 

U03-40 

169* 


J 1 = J 1 + 1 

U03M 1 

1 7u* 


8 J 1 =B J I ♦ I . 

003-42 

< 7 i * 


KTEK = Nfi-Jl-6 

CU3-43 

172* 


F 4 I J ) =SUMOE V 

003-4*4 

173* 


AH ( 1 , I J ) =A 1 

C03-45 

17-4* 


AR ( 2, I J ) =6 l 

003-46 

175* 


AK ( 3 , I J ) = A 2 

003-47 

176* 


AR < M , I J I =B2 

00350 

177* 


AR 4 5 , I J ) =A3 

0035 1 

178* 


AR ( 6 , I J ) =B3 

303 5 2 

17V* 


NKOT= I J 

303 5 3 

16c* 


IF 1kTER.EQ.0> 60 TO 305 

00355 

181* 


GO To 3300 

00356 

182* 

305 

CONTINUE 

00356 

183* 

C 
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OflASA 189* C SFtFcTTQN QF AhiD cnR-^SEtlbmj-N G — PAR-A MEIERS 


00356 

185* 

00357 

- 186* 

00362 

187* 

00369 

18 8* 

00365 

189* 

-0036®"“ 

19o* 


c 

DO 307 i» I ■ I d 

IF (F( 1 ) .LT.O. I GO To 306 

GO TO' 307 

3u6 F ( I > *-F ( I ) 

3xr7 CONTrtTtrE 


00370 

0037 1 

00372 
003 7 3 
0037** 
60375" 
00376 
OS **00 
00901 
00H02 
00*403 


19|* 
1 y 2* 
193 * 
19 * 4 * 
1 95 * 

- I fl l m 

170* 

197 * 
198 * 
199 * 
200 * 
2q 1 * 


JJl*| 

N K 0 Tr MX 0 T - 1 
J*1 

K *2 

TEMP ( I >*F ( 1 > 


IF (OF.GT«G») GO TO 3fl9 

K*K*1 — 

JJ1-K-1 

NKOT*NKOT-"t" 

IF (NKOT.EQ.fl) GO TO 3 I 0 


CO 90 5 2 02* GO TO 308 

CC*406 203* 309 J*JJl*l 

C0**07 20*4* TEMPI JT-FTK ) 

00*410 205* K *K ♦ 1 


00911 2Q6* JJl*K-l 

00*412 2C7* NKOT-NKOT-l 

00*4 1 3 20®* I F (^ OT *^e ri OT-1IO--T43- Jtfr 

00*4 1 5 2fl9* GO To 3i)8 

00*416 2 10* 3*G~ CONTINUE ~~ 

C0*4 1 7 2U* At«AR(l,J) 

00920 2 1 2*-- - B1«AR(2.J> 




00*42 1 213* A2*AR(3»J) 

0 0 *4 2 2 2 1 8* - - -8- 2»A R-t9r O ) 


00923 
009 29 

00925 

00926 
00930 

215* 
216* 
2 1 7* 
2 1 8* 
219* 



A3*AR (5, J) 

03*AH ( o i J ) 

JSEL3* J 

write (6-,-teo) - 

WRITE (6,101) 

®rO®32 



w •» i ft to»ss?; 

00939 

221* 


WRITE (6,120) 

009 36 

222* 

— 

DO “3 IT— T* 1 . 1 J — 

009*4 1 

223* 

311 

WRITE (6,121) ( I ,F < I > , ( AR ( J » I ) , J*l ,6 > > 

0095 2 

229* 

- ' 

Wft ITE T6jT0H - 

00959 

225* 


WRITE (6,122) 




* R I TE ( 6 i 123) JSC L 3 

00961 

227* 

90flC 

CONTI NUE 

00962 

228# 

‘ ' 

If ( ITERM3.EQ.C) GO To SO&C 

00962 

229* 

c 


00962 

230* 

- C INlTTAt ESTIMATES FOR T*0 EXPONENTIALS 

00962 

231* 

C 

STRAIGHT LINE FIT TO U 0G Y Vs. TIME 

"00962“ 

^32* 

V ' 


00969 

233* 


i «J*fc 

00966 

239*“ 


-J t » 2 

00966 

235* 


B J 1*3. 

00967 

236* 

901 

“ 

00970 

237* 


SM1*C. 

009 7 T 

233* 


SM2*0 • 

00972 

239* 


SM3*0« 
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G£j873 
fefe *4 7 8 
u£87S 
00500 
-66-5-0-1- 
66502 
00503 

00505 
005U6 

00506 
00506 
00506 
66507 
Jit 5 1 0 
035 1 t 
0051 i 
605 l 1 
00 5 i 1 
6C512 
005 1 3 
0051*1 
005 l 5 
&&5 1 6 
00517 
0&52C 
00523 
0052H 

00525 

00526 
w 052 7 
66536 
00531 
00533 
0053*1 
v 05 3*1 
0053*1 
0553*1 
005 35 
00536 

00536 
Cf 0 5 3 6 
665 3 6 

00537 
005*12 
055*13 
005*15 
005*16 

0555 1 

00552 

00553 
©»5S5 

00556 

00557 
feC 5 60 
6956 1 
00562 


2*10*. 
2*1 1 « 
2*12* 
213* 


581*3*- - - - 
J J * N N - J 1 
DO 1 * J J « N N 

SM 1 * 1 ( I ) • * 2 ♦ S M 1 

2 M *1 • 


S M 2 ° D"L 0 G ( Y ( I ) ) *SM2 

2*15* 


SH3 = T ( 1 ) ♦ SM 3 

2*46* 

*li»o 

SH*1*T<I)*0L0g(YU)>*SmH 

2*17* 


CAA2*BJl*SMl-SMi**2 

2*46* 


CA2«CSMi*SM2-SM3*SM‘ , )/CAA2 

2*49* 

C 

- - - 

2 SO* 

C ESTIMATES FOR A 2 i 8 2 

251* 

c 


252* 


B2 = -< 18J1*SM*4-SH3*SM2)/CAa2> 

253* 


A2*£xP ( C A 2 ) 

25*4* 


J2»Jl*3 

255* 

c 


256* 

C STRAIGHT line fit to log residual vs. 

257* 

c 


258* 


jn=nn- j 1 - » 

259* 


832" JN 

260* 


SM1«C. 

26 1* 


S M 2 ■ 0 • 

262* 


S M 3 * C* * 

263* 


Sfri*»*C. 

26*1* 


DO 8c 2 1=1, JN 

26b* 


X 2 ( l ) =-8 2* T t I ) 

266* 


KES(l)aY!I)-A2*DEXP<X3U)) 

267* 


RES ( 1 ) * D AB S t KtSI 1 ) ) 

268* 


SM 1 *T ( I ) **2*SM I 

2 69* 


SM2*DL0G ( RES t i ) ) +SM2 

271-* 


SM3 = T 1 I ) *SM3 

27 1 • 

*102 

SM9*T(I)*DL0g1RES(1I)*SM8 

27 2* 


CAA1*BJ2*SM1-SM3**2 

27 3* 


CA1*(SM1*SM2-SM3*SM8)/CAA1 

27*1* 

c 


275* 

C LSTI 

mates FOR A 1 , b 1 

276* 

c 


277* 


Bl«-UBJ2*SM*4-SM3*SM2)/CAAl > 

278* 


A 1 =EXP I C A 1 ) 

279* 

c 


2au* 

c weighted sum of squares or DEVIATIONS 

28 1 • 

c 


282* 


DO 833 1 = 1, Nil 

2 83* 


X 1 < 1 ) =-B I • T t I > 

2e*l* 

8*3 

X 2 ( I ) =-B2*T( 1 ) 

285* 


sumdevo. 

286* 


DO 8^.*4 1 = 1, Nn 

267* 


B l I , l ) * DE X P ( X l t I) ) 

288* 


81 I , 3 1 *DE XP 1 x 2 ( I ) ) 

289* 

Mel 

5UHDEV=<Y«J>-Ml*b(l*l)-A2*B t l»3>) 

29u* 


Jl=Jl+l 

291 * 


b J 1 = b J l ♦ 1_ . 

292* 


KTER=NN-Jl-3 

293* 


F l I 3 ) = S U 11 0 E V 

29*1* 


AR ( 1 , I J ) »A 1 

295* 


6 R ( 2', I J ) = B 1 


T 1 ME 
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- 0©54^3 2.9 6_* AK(3,1 Jl*A2 


00568 297* 

03565 2 96* 

00566 299* 

60570 300* 

03571 3q1« 

0*|5 7l -302* 

0357 1 3q 3 • 

00571 3os* 

00572 305* 

00575 306* 

00577 3y 7 * 


AK ( 8, I J) =B2 
NKOT * I J 

IF (KTER.EQ.OJ GO TO 435 
GO TO SOI 
Rub CONTINUE 

c selection of smallest f ano corresponding Parameters 

c 

DO Sy7 I * 1 » I J 

IF (Fill . L T • u • ) GO To 806 

GO To SO 7 


00601 


30«» 

309* 


St rto 
Hoi 


00603 3iU* 

0060S 3 1 i • 

00605 312 * 

00606 313 * 

CC607 3 1 8 • 

00610 3 1 S • Sob 


CONTINUE 

JJl«l 

NK0T*NK0T-1 
J=i 
K * 2 

DF«TEMP( Jl-F (K 1 


C061 1 316* 

0C6 1 3 317* 

006 1 S 3 1 8 • 

0061 5 319* 

"Ot)4t6 — 320* 

00620 321* 

00621 32 2* 

00622 323* 

00623 328* 

006 2 S 325* 

6062 S 326*-- 

£0626 327* 

006 30 326* 

00631 329* 

00632 330* 

00633 331* 

00638 332* 


IF ( OF . GT". C • 1 GO TO Sq 9 
K«=K* 1 
JJl*K-l ■— 

NKOT-NKOT-l 


S 10 


GO TO 808 
Suy J*JJ| ♦ 1 

TEMP ( J ) *F ( K 1 
K = K + 1 — --- - 
JJI «K- 1 

NKO-T*N KOT -1 

IF (NKOT.EQ.Q) GO to SID 
- GO To SO 8 

8l« continue 

A l *AR ( 1,31- 
6 1 *AR ( 2 i J 1 

A 2»AR t 3 i Ji - 



VU635 333* 

00636 33S* 

00637 335* 

G96S1 336* 

00683 337* 

00685 338 » 

60687 339* 

00652 3S0* 

00663 3 S 1 * 

00665 - 3S2* 
00667 3S3* 

— 006 7 2 3S-S* 


B2*AR(8,J) 

JSEL2* J 
WRITE (6,100) 

WRITE (6,101) 

WRITE ( 6 , l 2 S ) 

WRI TE (6,12 5 ) 

DO Si 1 1*1 , I J 

S l 1- WR LTE— ( 6-, 1 26 )— ( I ,f ( l ) , ( AR I J , I ) ,0*1 ,S ) ) 

WRITE (6,101) 

- „VRITE-(-6tI22) 

WRITE (6,123) JSEL2 

5u03 C ONTINU E 


0Q672 

3S5* 

03672 

3S6* 

33672 

3 S 7 • 

60672 

3S8* 

006 7 3 

3s9* 

00675 — 

— 3Su* 

00675 

351* 


C 

C STATEMENTS NUWbEREO 5^0"599 REFER TO PROGRAM FOR THREE 
C EXPONENTIALS 

C - ... 

IF (ITERMS.Eq. 01 GO To 600 
■560 CONTINUE — - - 

c 
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££6 7 S 
£3675 
£0675 
£3 6 7 6 
£0-7.0 0_ 
„j7 j2 
£37J1 

££706 

£j7i;j 

£0725 

££7 2 I 

£0722 

53723 

&072H 

WO 7 25 

£3726 

037 2 7 

0073u 

0373 1 

£073 3 

C373S 

05737 

£.£712 

0(J713 

C3711 

007 '15 

03 7*16 

ti£ 7 1 7 

£3755 

£575 1 

i!5 7S2 

£5753 

0375*1 

£3755 

£3/56 

£3757 

00765 

OJ766 

3577 1 

£3772 

03773 

£0771 

0377*1 

£377*1 

00771 

£3775 

01330 

£1550 

5 1 OGU 

£1530 

£1032 

£153 3 

01 30*1 

£ 1556 

£1311 

31311 


3s£* 

353* 
351* 
355* 
356* 
357* 
358* 
359* 
3 6 u* 
361* 
362* 
36 3* 
361* 

36 5* 
366* 
367* 
368* 
369* 
373* 

37 1* 
372* 
373* 
371* 
3 7 5* 
376* 
377* 
378* 
379* 
385* 

38 1 * 
382* 
383* 
381* 
385* 

38 6* 
387* 
388* 
389* 
395* 

39 1 » 
392* 
393* 
391* 
395* 
396* 
397* 
398* 
399* 
100 * 
lul • 
102* 
103* 
1£1* 
105* 
106* 
107* 


C THU El EXPONENTIALS 
C [1ERM1=l 

c 

WRITE (6| luU) 

-.1*1 T E ( 6, 1.0.2 ) 

N KITE (6,101) 

WRITE ( 6 , I y 3 ) 

(-•RITE (6,101) 

.(RITE (6,105) a 1 »8t »A2,B2 ,a 3,B3 
KCYCLE=0 
A i l = A I 
811*81 
A 2 2* A 2 
822=82 
A 3 3* A 3 
833=83 
8 J = N N 
50 MS W= J • 

WRITE 16,112) 

WRITE ( 6 , 1 jl ) 

WRITE (6,110) 

0 0 5 00 1 1*1 .Nil 

x n n«-a i i*t { 1 1 

X 2 ( 1 ) * -8 2 2 * T ( ( ) 

X 3 I 1 )s-833*T ( 1 ) 

8(1,1) =0EXP ( X 1 ( I 1 ) 

B ( I , 31 *OEXP ( X2( I ) ) 

8 ( I ,S1*0EXP(X3( T) 1 

Y C I I ) = A 1 1 • B < I ,1 )+A22*B( I ,3) + A3 3*8( I ,5) 

RES ( I ) *Y ( 1 I -YC ( I I 

R SS ( I ) *RES( 1 ) 

TSM )*T( II 

Y 5 ( I ) * Y ( I I 

Y C S ( I 1 ■ Y C ( I » 

SUMSlj = RE5< 1 )**2*SUMSy 
5gjl WRITE (6,128) Y ( 1) , YC( I ) ,RES( i ) 
SIGYlN*SOKT(SOHSQ/(BJ-6.)) 

WRITE (6,113) SI5Y1M 
X X = N N 

CALL 0 U I X 3 L l-I ,TL,TK,YB,YT.13,BCDT,BCDY,KK*TS r YS) 
CALL QUIK3L ( U, Tl - , Y8 , YT , 35 ,BCDT ,BCDY ,-KX , T5 , YC5 ) 
CALL QU1X3L<-I,TL,TR,Y8I,YT1,35,8CDT,BCDK,-KX,TS,RSS) 
C 

C A H A T K I X 

c 

00 5 y | 1 = 1 , NN 

Sul -7 l I ) * I . / V A R Y 

c 

C 8 MATRIX 
C 

5ol 0 CONTINUE 

XCYCLE-KCYCLE*! 

WRITE (6,1011 
WRITE (6,106) KCYCLt 
DO 5 52 1*1 , N (I 
X l ( I )=-B U*T { 1 ) 
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GlOlS HilO* X2( 1 ) s-B22*r (i ) 

Olul6 9Q9* 5i>2 X J < 1 > =-033*T ( I) 

21u20 H|u« WRITE ( 6 , 1 3 8 j 

01 02 2 00 5-J3 1 3 i ,HN 

ii 1 32 5 ‘UZ* 3 < I i t J *0EXP < X I ( I J ) 

0102 6 9 i 3 • B ( I ,2)*-All*T( l )*B(1 t j ) 

313 27 919* U( I ,3>=0£XP<X2( l )» 

0 103U His* 3 ( I ,9) *-A22* r ( I ) «& ( I , 3 > 

01031 **16* »( ! ,S)=DeXP(X3< 1 ) ) 

0 1<232 *♦ 1 7 • B ( I , 6 ) =-A33*T (1 ) *3 ( I , S ) 

01033 H10* Su3 WRITE (6,123) < B ( I , J B ) , J B = 1 , 6 > 

01033 919* C 

01033 92u* C TRANSPOSE OF B MATRIX 

**1033 921* C 

illO'42 922* DO 5 U 9 K* I , 6 

W 1 o 9 5 9 2 3 * 0 0 S J 9 1*1 » N N 

01050 929* SU9 BT (K , i ) *B ( I ,K ) 

01050 925* C 

OICSO *2 6* C N MATRIX 

01050 127* C 

31053 428* DO 5.JS 1*1, Nn 

G 1 05 6 929* 5 *,5 I3N ( l , 1 ) * Y l I > -A 1 1 *B ( I , | ) -A22*B ( I , 3 ) -A33*B ( I , 5 ) 

01056 930* C 

0 1 056 931* C B-TKANSPOSE*W*N 

01056 932# C 

OlLitoU 933* DO 506 1*1 , NN 

31063 939* 566 D 1 ( 1 , 1 ) =0 • u D 0 

31065 93t,» DO 5tj 7 N * 1 , N i ( 

01070 936* 5u7 OHN,| )*tf(N»*BN(N,l> 

31072 937* DO 508 1=1,6 

G 1 075 938* 5ub 02 < I , 1 ) *:j. UDq 

01U77 939* DO 5o9 M * 1 , 6 

01102 980* DO 5u9 N*1 ,Nn 

01105 981* 5y9 D2(M, 1 )*BT(M,N) *D| (N, j )+D2(M, I ) 

01105 992* C 

w 1 I u5 9.43* c ti-TRAN5P0SE*W*b 

Oil u5 899* C 

OHIO 995* DO 5 1 -3 1*1, Nn 

01 1 13 9 9 e> • DO 5 jq J=1 ,6 

01U6 987* 5 1 y D| ( 1 , J)*0*u D U 

01 121 998* DO S)t P*1 ,6 

01 129 999* DO 5)1 N * I , N N 

31127 9 SO* 51 1 Dl (N,P)*WtN)*ri(N,P> 

01132 951* DO 5)2 1*1,6 

01 135 952* DO 512 J=l ,6 

01 190 953* 512 D3( I , J)*G.OOo 

01183 959* DO 513 M*l,6 

01 196 955* DO 513 P*1 ,6 

31151 956* DO 513 N*l,rtN 

01159 957* 513 D3 ( M , P ) *BT 1 M , w ) *0 1 ( N , p ) *Q 3 ( N , P ) 

£1160 958* WRITE (6,199) 

01162 959* DO 5133 1*1,6 

C 1 l 65 960* Si3y WRITE ( 6 , 1 2 8 ) ( D 3 ( 1 • I X ) , I X * l , 6 ) 

£1165 961* C 

£1165 9fc2* C INVERSE OF B - T R A N S P 0 SE * « • B 

01165 963* C 
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©1174 
C 1 1 7 5 
w l 200 
V# 1 2 U 3 
O' l 2 0 4 

464# . 

465# 
466* 
467* 

468 • 


12*0 , 

DO 5 | 4 1 * 1,6 
00 514 J*1 ,6 
12 * 12*1 
D ( 1 2J ■ D 3 LJ , I ) 


4fc9* 

514 

DO 1 I 2 ) *D 3 ( J • 1 > 

61 210 

470* 



N* 6 — ' '■ * ’ ~~ 

a» 2 i i 

471* 


M = 0 

©1212 

472* 


C ALL I N V R T ( D , N i M , DETe« ) 

V i 2 1 2 

473* 

C 


©1212 

474* 

C 1 N VERSE MATRIX ------ 

& 1 2 1 2 

475* 

C 


01213 

476* 


Kl« 1 

«. 01214 

477* 


K 2“ 6 

01215 

478* 


K 3 * 6 

01216 

479* 


K 4 * 6 

01217 

48o* 


WRITE (6,137) 

01221 

481* 

Si Hu 

WRITE (6,129) ( 0 ( I ) , 1 =K l , K 2 * 

~ 0t227- 

- 482* 


K 1 * K 2 * 1 

" 01230 

483* 


K2=K2*K3 

01231 

464* 


K 4 * K 4 - I 

01232 

485* 


If (K4.EW.0l 00 To 5 1 4 1 

01234 

486* 


GO TO 5140 

01235 

487* 

5141 

CONTINUE 

"01235 

488* 

c 

. . 

01235 

489* 

C 1NVERSE*0R1GINAL MATRIX 

01235 

490* 

C 


01236 

491* 


12 *j 

01237 

492* 


00 5(5 1 * 1,6 

01242 

493* 


DO 5 | 5 0*1 ,6 

01245 

499* — 



12 * 1 2 * I 

01246 

495* 

515 

0 3 3 ( J , I ) *0 ( I 2 ) 

01251 

496* 


DO 5 1 50 1 * 1*6 

01254 

497* 


DO 5 1 5Q J*1 ,6 

01257 

498* 

5150 

D 1 DEN ( I , J)=U.cD© 

01262 

499* 


DO 5151 M*1 ,6 

01265 

500* 


DO 5151 P* l . 6 

01 270 

So 1 • 


DO 5151 N* 1 , 6 

01273 

502* 

5151 

DlDENCM,P)*D33(M,NJ*D3<N,P)*DIDEN(M,f») 

01277 

503* 


1 n*6 

01300 

5C4* 


1*1 

01301 

SC5* 


WRITE (6,127) 

01303 

53** — 

5152 " 

KR t rr ( 6 , |28r TDIOENl I , J ) , J* 1 , 6 ) 

01311 

507* 


I W* I ft- I 

01312 

5Q8* 


IF 1 I (V • E W * 0 1 GO TO 5js3 

01314 

509* 


1*1 + 1 

01315 

510* 


GO TO 5152 

01316 

511* 

5153 

continue 

01317 

5 1 2* 

— 

(WRITE ( 6 , 10»» - -- 

01321 

5 13* 


WRITE (6,129) 

01323 

514* 


WRITE (6,128) DETER 

01323 

5 IS* 

C 


C 1 323 

b l 6 * 

C DELTA MATRIX 

01323 

517* 

C 


01326 

5 I 8 * 


DO 516 1*1,6 

01331 

5 l 9 • 

516 

DEL ( I , 1 1 *C.UD*> 
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V » J J J 

6*1336 



S2i* 

01341 

522* 

01341 

523* 

01341 

524* 

01341 

525* 

— 

526* 

013 4 5 

52 7* 

0 1 3 4 6 

528* 

01347 

S 2 9 • 

0 350 

530* 

01351 

531* 

Cl 3ST 

- S'32~* 

C 1 35 I 

533* 

C 1 35 1 

534* 

S 1 352 

535* 

01355 

536* 

01356 

537* 

CT357 

5 38* 

01361 

539* 

01364 

540* 

01365 

54 1 • 

01366 

542* 

01370 

543* 

-CT372- - 

544* 

01374 

545* 

01376 

546* 

01 400 

547* 

i, 1 403 

548* 

01404 

549* 

0^05 

556* 

01406 

551* 

01407 

552. 

01407 

553* 

01407 

554* 

01407 

555* 

0*1 5 

566* 

01416 

557* 

01417 

558* 

01422 

55 9* 

01424 

560* 

01424 

56 1 • 

-0142-4 - 

--66-Z-* 

01424 

563* 

01425 

564* 

01426 

565* 

01431 

S66* 

014 32 

56 7* 

01 432 

568^” 

C 1 4 32 

569* 

C 1 4 3 2 

570* 

01434 

571* 

£1435 

572* 

01436 

573* 

- C l 4 37 

5 7 4* 

01440 

575* 


DO 517 N*1 ,6 

517 DEL<M,1>»D33<M,N1*D2(n,I)*DEL<M,1) 

C 

C IMPROVED PARAMETER ESTIMATES 
C 

*t-t»|ttT*OeL« | f| ) " — •• “ ' ' 

8 1 1 »B 1 | + DEL • 2 i 1 ) 

A22»A22*OEL (3.1) 

822“B22*DEL ( Mil) 

A33*A33*DEL<$, 1 ) 

B33*B33*DEL ( 6.1) 

C RESIDUALS USING IMPROVED PARAMETER ESTIMATES 
C 

DO 519 1*1 ,NN 
II I I l**B11*T III 
X2 ( t ) «-822»T | 1 ) 

5T9 X3Tn*-0 33*T rn ■“ ' 

DO 520 I'i.NN 
B ( I • I ) »D E X P ( X 1 ( I ) ) 

B( I ,3)*DEXP(X2( I ) ) 

S2l B ( I 1 5 ) * D E X P I X 3 ( I ) ) 

WRITE 16,100) 

: — ion Te-rerrrr* r~ — 

WRITE (6,101) 

WRITE 16,140 
DO 521 I > I , NN 

- YC< I )=A1 1 *8 ( I , | ) ♦ A 22 *B ( I ,3)+A33*B( I .5) 

RES ( I 1 - Y ( I ) -YC I 1 ) 

*ss-<-n •ittsi-H- 

yes < 1 1 ■ y c 1 1 ) 

52 1 write (6,128) Y ( i ) , Y C ( I I , RE S ( I ) 

C 

C STANDARD DEVIATION OF RESPONSE VARIABLE USING INITIAL ESTIMATES 

c 

SUM1« 0 , 

DO 522 I *1 ,NN 

522 SUM1«BN( I ,1 )««2*SUMI 
S1SY 1 «SQRT (Sum / (BJ-6, ) ) 

C 

■-C— ST^ANO"ARfr OEV I AT 1 Oh OF R5SP0N5E VARIABLE US ) NG “TMPStrVeir EST t M A TE S 
C 

SUM2«C* 

DO 523 I«l ,NN 

5 1GY2«S0RT ( SUM2/ ( B J-6. ) ) . . 

523 SUM2 = RES( I )**2+SUM2 

C PARAMETER STANDARD DEVIATIONS 
C 

SGA1=SIGY2/SQRT(D3(I,1>) 

5GBI*SIGy 2/S QRT (03(2,2) ) 

SGA2=S!GY2/5QRT(03(3,3)l 

SG-8-2T=SrG Y2AStfTrrTD3 fM-rM ) ) — — 

SGA3«SIGY2/SQRT (03(5,5) ) 
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61 ** * 
is l * * 2 
l> 1 * * * 
01**6 
C.1_*5JL 
«1 *56 
01*60 
91 H 7 C 
£1*72 
W l H 7 M 
oisoo 
5 1 SO 2 
S1S0 1 * 
£150* 
Si l 50* 
w * SO *< 
0151* 
(si s |5 
01 S I 6 
C'1S»7 
01521 

LIS22 
C 1 523 
CIS2S 

015 30 
01S33 
01S3S 
01637 
tl 5 < t 2 
CIS** 
C1S*«> 
015*7 
C 1 560 
C i 56 1 

01562 
01553 
C 1 563 

01563 
01653 
£1563 
01565 
01567 
01567 
01567 
01557 
is 1 56 7 
OlS&C 
01662 
£ l 66 * 
£1566 
01570 
01572 
C160C 
01601 
51602 
01603 


576* 
577. 
57 8* 
679* 
580. 
581* 
582* 
583* 
58*. 
685. 
586* 
587* 
588. 
669. 
590* 
59 1 • 
592* 
593* 
59*. 
595* 
596. 
697* 
598* 
599. 
600 * 
60s* 
6 £ 2 * 
603. 
60** 
60S* 
6C*. 
607* 
6 C 8 * 
6 0 9 • 
6 1 0* 
6 1 1 . 
6 l 2» 
6 1 3* 
6 1*. 
615* 
6 16 * 
6 l 7 * 

6 l 8. 
6 1 9. 
6 20-. 
621* 
6 22 * 
623* 
62 *. 
625* 
626. 
627. 
628* 
629. 
6 3C * 
63 I * 


S6B3*SI6Y2/St,jKT ( 03 ( 6 ,6 ) ) 

WRITE. (6,101) 

WRITE (6,lt‘7) 

WRITE (6,106) (DEH l . I ) . 1 = 1 .6 I 

SV_R IT E__ ( 6,1 Cl) 

WRITE (6,109) 

WRITE (6,t(j5| A I 1 ,8) 1 , A 2 2 ,B22,A33,B33 
WRITE ( 6 , 1 u I ) 

WRITE (6,110) 

WRITE (6,111! S I GY l . S I GY2 
WRITE (6,101) 

‘ WR HE ( 6 , II 2) 

WRITE (6,1 13) SGA|,SGB1»SGA2,SGB2,SGA3,S0B3 
C 

C iteration LOGIC CYCLE 

c 

T0LI=SIGY2»*2-SIGYI*»2 

OTOL*TOLER“TCiLl 

DTOL-=ABS(OTOl j 

IF ( ( OTOL-.OOu&l ) • GT .£,. ) GO TO 52* 

GO TO 525 
52* T0LER«T0L1 

WRITE ( 6 , t C 1 ) 

WRITE (6,11*) TOLER 
WRITE (6,1*1! KCYCLE 
525 WRITE ( 6 , 1 0 1 i 
WRITE (6,100) 

WRITE (6,115) DTOL 
WRITE (6,116) 

WRITE (6, II 7) 

KK=WN - - — 

CALL QUIK3L (- l , Tl ,TR , YB , YT . *3 ,BCDT ,BCDY ,KK , T5 , YS ) 

CALL 0 U I K 3 L ( U , TL , TR , YB , YT , 35 , BCDT ,BCD Y ,-KK , TS , YCS ) 
CALL QUIK3L I- 1 ,, TL , TR , YB I , YT l , 35 ,BCDT ,BCDR ,-KK ,TS ,RS5 ) 
NCASESbNCASES- 1 
IF (NCASES.EW.C) GO to 800 
C 

C STATEMENTS NUMBERED 60c" 699 REFER TO PROGRAM FOR TWO 
C EXPONENTIALS 

c 

IF ( I TERM3 . E(j .C > GO To 200 
6uC CONTINUE 
C 

C TWO EXPONENTIALS 
C I V E R M 3 = I 

C 

WRITE (6, ICC) 

WRITE (6,102! 

WRITE (6,101) 

WRITE (6,130) 

WRITE 16,10*) 

WRITE (6,131) AJ.BI.A2, 82 
KCYCLE.O 
A I I -A I 

811=81 ' 

A 2 2 = A 2 
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01604 
U 1 605 
u 1 606 
il6w«7 
w 1 6 11 
u 1 6 1 3 
a I 6 I 5 
0162C 
0 1 62 I 
01622 
£1623 

a 1 6 2*t 

{,'1625 

01626 

01627 

S163C 

01631 

81632 

81633 
01641 
£16 4 2 
81645 

01646 

01647 
£ 1 6 5 1 
01650 
01650 
Cl 650 
u 1 65 1 
01654 
01654 
01654 
01654 

01656 

01657 
016 60 
01662 
01665 

0 1 670 
01671 
81673 
016 7 5 
£1700 

01 701 

01702 

01703 
SI 704 
0170* 

0 l 7 C 4 

01 704 
£ 1 7 1 3 
01 7 1 6 
01721 
£ 1 7 2 1 
01721 
01721 


632* 


622=82 

633* 


B J * N N 

634* 


SUMS(,' = £. 

635* 


WRITE (6,142) 

636* 


WRITE ( 6 , 1 0. 1 ) 

637* 


i>, R I T E ( 6 , 1 4 1 ) 

63b* 


DO 6 COI 1*1. NT. 

63 V* 


X 1 ( I >=-B | 1*T ( 1 ) 

6 4c* 


X 2 ( I ) = -8 22 • T ( 1 ) 

64 i • 


6(1,1) = DEXP ! X | ( I ) ) 

642* 


6 < 1 »3)«DEXp(X2( I ) ) 

643* 


YC< 1 )=AI 1*6( I ,1 )+A22*b< 1,3) 

644* 


RES ( I ) *Y ( I) -YC ( I) 

645* 


RSS ( 1 ) =RES ( I ) 

646* 


TSC 1 ) = T< I ) 

647* 


YS< I > = Y ( I ) 

64b* 


YCSm«rYC(M 

649* 


SUMS(,'*RES( 1 )**2*5UMSu 

650 * 

6001 

WRITE (6,128) YU ) ,YC( I ) ,RES 

651 • 


SI6YIN*SQRT(S0MSQ/(Bj-4.) ) 

652* 


WRITE (6,143) Si GY IN 

653* 


KK«NN 

654* 


CALL QUIK3L (-1 ,TL,TR,YB,YT . 

655* 


CALL QUIK3L (0,TL,TK,YB,YT ,3 

656* 


CALL QUI K3L ( - 1 , TL , TR , YB 1 , YT 1 

657* 

c 


6 5 8 * 

C iv 

MATRIX 

659* 

C 


660* 


DO 6c 1 I = l , NN 

66 1 • 

6ul 

W( I )*( ,/VARY 

662* 

C 


663* 

C U 

MATRIX 

664* 

c 


665* 

6c 1 v. 

CONTINUE 

666* 


KCYCLE*KCYCLE* 1 

667* 


W KITE (6,101) 

66b* 


WRITE ( 6 , 106 ) KC YCLC 

669* 


DO 602 I * 1 , Nn 

67c* 


XIII) =-B 1 1 *T ( ) ) 

671* 

6u2 

X 2 ( I 1 * -B 2 2* T ( 1 ) 

672* 


WRITE (6,13b) 

673* 


DO 6i,3 1 = 1. NN 

674* 


B ( I . 1 ) =DEXP ( X 1 ( I ) ) 

675* 


fllI.2l=-A||*Tlll*B|l,|] 

676* 


8(1,3) = 0 E X P t X2 ( I 1 ) 

677* 


B l I .4) =-A22*T ( I 1 *b ( I , 3 » 

678* 

6 v3 

WHITE (6,12b) l B ( 1 , JB ) , JB= ) , 

679* 

C 


6 8 0 * 

C TRANSPOSE OF B MATRIX 

68 l * 

C 


682* 


DO 604 K = 1 , 4 

683* 


DO 6 o 4 1*1 , N N 

684* 

604 

8T (K , I) = B ( I ,.K ) 

665* 

C 


686* 

C N 

MATRIX 

687* 

C 



■KK .TS,YCS> 


4 ) 
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Cl 72*4 
01727 
01727 
01727 

686* 
689* 
690* 
6 V 1 • 

6i.5 

C 

C B- 

DQ 6©S -I-S'-i .UN - 

BN<I,l)=r<l>-All*B(l,i>-A22*l 

TRANSP0SE*W*N 

01727 

692* 

C 


017 31 

693* 


DO 6c 6 I«1,NN 

Cl 73*4 

6 9 9* 

606 

oit r, i )*c*cdo 

01 736 

695* 


00 607 N* 1 , NN 

Cl 7M1 

6 96* 

6u7 

DI<N,l)«W4M)*BN(N,l) 

0 1 7 *» 3 

697* 


00 6 08 1*1, *4 

01 7 96 

698* 

60b 

02(1, 1 >«O.CDo 

01750 

699* 


DO 609 M aa 1 , *4 

31 753 

700* 


' DO 6C.9 N«I ,NN 

01756 

701* 

609 

U 2 ( M , | ) * B T ( M , N ) * 0 1 (N, , > ♦ D 2 ( H 

Cl 756 

7 C2 • 

C 


0175 6 

7 0 3 * 

C b- 

transpose *w*b 

01756 

7 0*4 * 

C 


017 6 1 

7 0S* 


U 0 6 10 I = I , N N 

Cl 76*4 

7 06 • 


DO 610 J=l ,*4 

01767 

7j7* 

610 

D 1 ( I , J ) «0*cDy 

01772 

708* 


DO 611 P*l ,*4 

01775 

709* 


0 0 6 J | N * | , N n 

w 2 000 

7 10* 

6» 1 

01 (N,P)=A(Nl*B( N,P) 

02003 

711* 


DO 6 J 2 I = 1 t*4 

32006 

7 12* 


DO 612 J * 1 , *4 

0 2C 1 1 

713* 

612 

D 3 ( I i J ) x v»u^O 

3201*4 

7 l H* 


DO 613 Maaj ,<4 

02017 

715* 


DO 6 | 3 P*i , *4 

02022 

716* 


DO 6 | 3 N-l.NN 

©2025 

717* 

6 1 3 

D3(M,P)»BT<M,|m)*01 ( N » P ) * D 3 ( M 

S2031 

718* 


WRITE T6, 1 *»*♦ r - 

02033 

7 19* 


DO 6J30 1 B l a *4 

©20*36 

720* 

6 130 

WRITE (6,128) (03(1, IX) • I X ** 1 

02036 

721* 

C 


32236 

722* 

C INVERSE OF B-TR ANSPOSE*W*b 

02036 

723* 

C 


020*45 

. 72*4* 


12*1 

020*46 

725* 


DO 6|9 1*1 , *4 

02051 

726* 


DO 6 | H J* I , *4 

0205*4 

727* 


D 4 ! 2 ) *03 ( j , I ) 

02055 

728* 

6 1 <4 

12*12*1 

0 2 0 6 C 

729* 


N * *4 

C 2 u 6 1 

730* 


m=o 

02062 

731* 


call INVR1 ( 0 , N , M , D t- T E R ) 

02062 

732* 

C 


62062 

733* 

C INVERSE MATRIX 

32062 

7'3*4* 

c 


32063 

735* 


Kl*l 

0206*4 

736* 

- - * 

K 2**4 ' 

62065 

737* 


K 3* *4 

02066 

738* 


K *4 * *4 

62067 

739* 


WRITE (6,137) 

C2071 

7 M C • 

6 1 *4C 

WRITE (6,128) (D( I ) » 1=KI ,K2> 

02077 

7 S 1 * 


M**2+ 1 

02 ICO 

7*42* 


K2*K2*K3 

32101 

7 *4 3 • 


1 

3“ 

XL 

W 

4“ 

* 



uzmz 

744« 

• - 

. 4.E_..iK4.*£0*0) GO 10 6(41 

32104 

745* 


GO TO 6146 

02105 

7 46* 

614 1 

continue 

02106 

74 7* 


12*0 

62107 

748* 

' " 

DO 615 1*1,4 

62112 

749* 


■00 615 0*1.4 

62115 

7 50 r" 


12*12+1 

02116 

751* 

615 

033< J • 1 > *D ( I 2 > 

62116 

752* 

C 


62116 

753* 

C INVERSE*0RJGINAL MATRIX 

62116 

754* 

C 


62121 

755* 


00 615C 1*1*4 

62125 

756* 

. 

00 6150 J*l *4 - — ' 

62127 

757* 

6 1 So 

0 I OEN ( I , J 1 * y . w . D 0 

02132 

758* 


00 6151 M=1 ,4 

62135 

759* 


DO 6J5I P* 1 • 4 

C 2 1 4 G 

760* 


DO 6151 N*l,4 

02143 

76 1* 

6151 

01 DENIM, P)*D33(M,N)*03<N,P)+01DEN(M,P) 

62147 

76 Z* 


I W*4 

62150 

763* 


1*1 

62151 

764* 


WRITE (6,127) 

62153 

765* 

6*52 

/.'RITE (6,12®) (DJDENl I ,J) ,J“i ,4) 

62161 

766* 


IW* I W-| 

62162 

767* 


IF (IW.EQ.O) GO TO 6153 

02164 

768* 

' - 

I»l*l 

62165 

769* 


GO TO 6152 

62166 

7 7 0* 

6 153 

continue 

62167 

77 1* 


WRITE (6,101) 

62171 

772* 


WRITE (6,129) 

62173 

773* 


WRITE (6,128) DETER 

62173 

7 7 4 * — 

C 

- 

02173 

775* 

C DELTA MATRIX 

62173 

776* 

C 


62176 

777* 


00 616 1*1,4 

62 201 

778* 

6 16 

DEL < I , 1 > *O.ODo 

62203 

779* 


DC 617 Mb l , 4 

62206 

7 60* 


00 617 N* 1 ,4 

6221 1 

78 i * 

6 1 7 

DEL<M,|)“D33(II,N)*D2(N,I)+DEu(M,1) 

6221 1 

782* 

C 


6221 1 

783* 

C IMPROVED PARAMETER ESTIMATES 

62211 

784* 

C 


£22 | 4 

785* 


A l 1 *A I l +DEL (1,1) 

02215 

786* 


8 1 1 *B 1 I +OEL (2,1) 

62216 

787* 


A22*A22+0EL (3,1) 

62217 

788* 


B22*B22+DEL ( 4 , 1 ) 

62217 

789* 

c 


62217 

790* 

C RESIDUALS using improved PARAMETER ESTIMATES 

62217 

79 1 * 

c 


62220 

792* 


D CT 6 | 9 ' I * 1 , N ft 

62223 

793* 


X 1 ( I ) *-B I i *T ( I ) 

62224 

794* 

619 

X 2 t I ) xs-B22*T ( I ) 

62226 

795* 


00 623 r=l.NN 

6223 1 

796* 


B ( I , 1 ) *OEXP ( X 1 ( I ) ) 

62232 

797* 

6 2 U 

8(1,3) *0£XP ( X 2 ( I ) ) 

62234 

798* 


WRITE ( 6 , 1 oc ) 

02236 

799* 


WRITE (6,139) 
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022*42 

022*4*1 

022*47 

022SC 


WRITE (6,1*40) 

00 621 1*1, NN 

YC(I)-A1I*B(1,1)*A22*b(I.3) 
REST I )«Y ( I )-YC ( I ) 


-02-25-1 

80S* 



RSS1 I ) • RE S < 1 ) 

C2"2S2 

806* 


TCS ( I ) * YC ( I 1 

02253 

8C7» 

62 | 

WRITE (6,128) Y( I ) ,TC ( 1 ) ,RES( I ) 

02253 

808* 

C 


022S3 

8C9* 

c standard deviation of response variable 

02 253 

8 1 u* 

c 


022 61 

8 | 1 • 


B J*NN 

02262 

812 * 


sum sir; 

C2263 

8(3* 


DO 622 1* 1 ,Nn 

02266 

81*4* 

622 

SUM|»BN( I , | ) •• 2+SUM | 

02270 

8)5* 


SIGYl*SfiRT ( SUM 1 / ( B J”*4 . ) ) 

02270 

8 16* 

C 


02270 

8 1 7 • 

C STANDARD DEVIATION OF response V a r,aBLE 

C227C 




02271 

819* 


S UM 2 = G • 

02272 

82C* 


DO 623 1*1 ,NN 

02275 

821* 

623 

SUM2*RES t I ) **2*SUH2 

02277 

822* 


S1GY2*S8RT(SUN2/(6J-5.)I 

02277 

823* 

c 


- 02277 

8 2 <4* C-pTRFTTET E R “ S TS74D A-R D~ D E V I A T I ON S 

02277 

825* 

c 


02300 

826* - 

_ 

SGA 1 «SI GY2/SQKT tD3( 1,1)1 

02301 

827* 


SGB1«S1GY2/SQRT(D3(2,2>> 

C23C2 

828* 


SGA2*SIGY2/SqRT ( D 3 ( 3 , 3 ) ) 

02303 

829* 


SGB2*SIGY2/SQRT(D3(*4,*4)) 

0230*4 — - 

83C* 

— 

- FRITE Tincn ' 

02306 

831* 


WRITE (6,107) 

02310 

832* 


WRITE (6,132) (DELI 1 , 1 > * 1*1 I *4) 

02316 

833* 


write ( 6,1 an 

02320 

83*4* 


WRITE (6,109) 

02322 

835* 


WRITE (6,133) A 1 1 ,B> 1 ,822,b22 

02330 

83 6* 


- Y.RITE t6i TCI 1 " 

02332 

837* 


WRITE (6,110) 

0233*4 

838* 


WRITE (6,111) SIGYI , S I G Y 2 

0 2 3*4 C 

839* 


WRITE (6,101) 

023*42 

8*40* 


WRITE (6,112) 

023*4*4 

8*41* 


WRITE (6,13**) SGA| ,SGBl ,SGA2,SG82 

023*4*4 ' 

8*12* 

c 


023*4*4 

8*43* 

c iteration logic cycle 

023*4*4 

8*4*)* 

c 

" - - 

02352 

8*45* 


T0L1»SIGT2«»2-SIGYI**2 

02353 

8*46* 


dtol»toler-tol I 

0235*1 

8*47* 


OTOL-ABS ( DTOl ) 

02355 - 

8*48* - - 


— t"F TTOTOL-vOCUal ) ,GT . a « ) GO TO -62*4 

02357 

8*49* 


GO TO 625 

02360 

850* 

62*4 

TOLER*TOL 1 

02361 

851* 


WRITE (6,101) 

02363 

852* 


WRITE (6,11*4) TOLER 

02366 

663* 


WRITE (6,1*41) KCYCLE 

02371 

-85*4* -- 

— 

WRITE (6YTCT0T ' 

02373 

855* 

625 

WRITE (6,101) 
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62375- 

6 56» . 

- _ 

- .WRITE (6,1154-0X06— - -- 

92-400 

857* 


WRITE (6,116) 

62-402 

856* 


WRITE (6,117) 

62-40-4 

859* 


kk»nn 

w2'fOS 

86c* 


CALL QUIK3L < - 1 , Tl . Tr , YB , r T . -4 3 , 8 CO T , BCD r , K K , TS , T S ) 

C2-4(,6 

861 ♦ 


CALL 0UIK3L I , TL , TK , Y6 , Y T . 35 , BCDT , BCD Y , -KR , TS , YC S ) 

£»2-4o7 

862* 


CALL CIUI K3L < -1 ,TL,TR,YBl.YTI,35,BCDT,BCDR,-KK,TS,RSS, 

62-4 1C 

863* 


WCAStS=NCAS£S- l 

62*41 | 

86-4* 


If INCASES. EC. t> GO To 

02-413 

865* 


GO TO 2w0 

02-413 

866* 

C 


C 2 -4 i 3 

867* 

c 

STATEMENTS NUMBERED 7C.O-799 REFER TO PROGRAM FOR ONE EXPONENTIAL 

02-4 13 

868* 

c 

I TEKM2= 1 

62-4 13 

869* 

c 


02-413 

876* 

c 


C2-41-4 

87 1 • 

7w u 

5M 1 *0 • jOO 

02-4 15 

872* 


S M 2 * 0 • C o 0 

62-416 

873* 


SM3*0 , (jOC 

02-4 17 

87-4* 


SM-4«0. JOG 

62-4 20 

875* 


B J 1 »NN 

62-421 

876* 


SUMSQ-C-. 

02-422 

877* 


DO 7(j | I * 1 , WN 

02-425 

878* 


SMI "T 11 ) *•* 2 ♦ S M 1 

62-426 

879* 


SM2*0L0G ( Y ( I ) ) ♦ SM2 

02*427 

880 * 

- 

SM3*T ( I ) +SM 3 

02-436 

88 1 • 

761 

SM -4 * T ( 1 1 *OLOG 4 Y ( 1 ) ) ♦ SW‘4 

02-432 

682* 


CAA 1 e BJl *SM 1-SM3**2 

62-433 

883* 


CAl*(SMl*SM2-SM3*SM-4)/CAAl 

02-43*4 

88*4* 


81 I « - < (Bji*SM‘4-SM3*Sm2>/CAA) ) 

62-4 35 

885* 


A 1 1 »DEXP ( CA 1 ) 

02-436 

886* 


WRITE (6,1 a 01 

02-4-46 

887* 


WRITE 4 6 , 1-15) 

02-4-42 

888* 


WRITE (6,1-46) 

6 2 -4 *4 *4 

889* 


WRITE (6,1-17) A 1 l ,B1 1 

02*4 50 

896* 


WRITE (6,139) 

62-452 

891* 


WRITE (6,101) 

02-45-4 

892* 


WRITE (6,1-40) 

02-456 

893* 


DO 7(,2 1*1, Nlv 

02-461 

89-4* 


X 1 ( 1 ) = -B 1 1 *T ( I ) 

02*462 

895* 


B ( I , 1 ) *DEXP 4 X 1 ( I ) ) 

02-4 6 3 

896* 


YC( I )=Al 1*B< 1 , 1 ) 

62-46-4 

897* 


RES ( 1 ) «Y ( I ) -YC ( I 1 

02-465 

898* 


KSS( I )*RES( I ) 

92-466 

899* 


TS( I )=T< II 

02-467 

906* 


YSI I )*Yl I ) 

62-470 

90 1 * 


YCS ( I I = YC ( I ) 

02-471 

9C2* 


SUM5Q «RE S ( I ) **2+SUMSc 

62-472 

903* 

7 4.2 

WRITE ( 6 , 128 ) Y ( I ) , TC | I ) ,R£S( I ) 

C2500 

90*4» 


SlGYIN*SQRT(5UMS(J/(8Jl-2. ) ) 

6250 1 

9£>5* 


WR I TE ( 6 , 1 -13 I S I G Y 1 N 

6250-4 

9 66 * 


KK«N(J 

025U5 

9(J7» 


CALL 0UIX3L < - l , T L , T R , Y B , Y T , H 3 , BCD T , B CD Y , K K , T S , Y S > 

62506 

9Q8 * 


CALL QUIK3L l w , TL , TR , YB , YT , 35 , B CD T , BC D Y , -K K , T S , Y C S ) 

62507 

9q9* 


CALL QUIK3L 4 - 1 , T L , T R , y B 1 , Y 1 1 , 35 , BC DT , BC DR , -K.K , T S , R Ss ) 

02510 

910* 


NCASES*NCASES - 1 

6251 1 

9 1 1 * 


IF INCASES. Eg. 01 GO To BOq 
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025 1 3 9 12* _ GCl-TO- 2tfQ 

&25IH 9j3* 80fc CALL ENDJOB 

025|5 9i*«* STOP 

it 2 5 } 6 9 l 5» END 


ENO OF COMPILATION: no DIAGNOSTICS, 
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epOK.iS Jmvhi , jwvnl 
HV 1 lb-135 32 !,«> 


SUbROUT I NL 1 N k it T El.IRf POInI 0vU36S 

storage; used: coDE(I) cc:R|3: datai„j c„oihr: blank cgmmoni 2) OuCOuu 


EXTERNAL REFERENCES (BLOCK. N a ml ) 

0003 tRacsf 
ObOR NtK«3» 


STORAGE ASSIGNMENT (8LOC*. TTPt, Rel* 1 ! Vfc LOCATION, NaHE) 


DuOi 

OOyu 1 b 11*6 

Ut'U 1 

UUvtRb 

1 l&L 

U&U 1 

Oouo27 »I7G 

00U» 

00003 A 

123b 

OOu 1 

•jQUuRZ 

1 3 »G 

OuUl 

Jwj 116 1 S7 q 

JLol 

J-..„2SS 

1 0 1 L 

w 0l 1 

G6ul*3 2C lb 

0001 

O0C203 

2?7g 

0001 

U0U2 2b 

220G 

0061 

3262*3 2 3 Rq 

u 66 1 

C 2-333 

235L 

-0-1 

GCo3(j7 2 hS(j 

O 06 J 

C003R3 

25CL 

0001 

J0C3RS 

25SL 

OOQl 

6Uji>71 7SL 

L'Cul 

G2Uj73 

76L 

jOuO 0 

Oi>-}CC2 AKAX 

CLOG i 

i COO 1 0 3 

1 

oocc 1 

1 000106 

1C 

bUO*J 1 

1 £uSl9* 1 NO 

SC UU 1 

[ G w J t; 1 

INDEX 

2030 1 

G60I10 I ND2 

*000 

0001 1 6 

INJPS 

0000 1 

oOOCOR 

IPI V 

DCu~ 1 

1 6Vc 1 L»S JR 

6L-w£ ! 

i wbU 1 62 

J 

jL-U 1 

OoOW)2 K 

lcl-c : 

1 uOOlOT 

L 

0000 1 

; -ooiet 

NN 


OGON D Qvzv'ju SIGN 


OOlQl 

1 * 



SUBROUTINE InvRT<a,N,m,0£T£K) 


06103 

2 * 



PARAMETER I D I M ■ 2u 


06 I U 3 

3 • 

C 


MATRIX INVERSION aN0 SIMULTANEOUS EQUATIONS SOLVER 


si J | 03 

R • 

c 


AMNPUT MATRIX FOR InvERISON OR AUGUrlENTED HaTHIX FOR S I mE . EWS. 

IGSCD'030 

66103 

S* 

c 


H»0«0EH of Coefficient matrix 

IGSOJORO 

0C103 

6 • 

c 


M B 0 FOR INVERlSON ONLY 


06103 

7* 

c 


M»NUMBER OF constant vectors 


CO 103 

6 • 

c 


DETEK.OETERM InaNT OF COEFFICIENT matrix 


GCIOR 

9. 



OOUB LE PRLClSlUN a < 2 1 . UETER , 5 / GN , AM A X 


601 OS 

10* 



0 1 MENS ION l P l V ( 1 0 l M 1 , l NO EM I 0 l M , 2 ) 


60106 

1 1 • 



DETER. | , COC 


03107 

12* 



SIGN. 1 ,000 


661 10 

13* 



DO 2 J J . 1 *14 


CGI 13 

I R • 


2S 

JPI VI J)*3 


031 15 

1 5 • 



NN.N.M 


00116 

16* 



DO 1 B 2 K«l ,N 


60121 

1 7 • 



AMAX»0. JOO 


60122 

IB* 


R3 

00 76 J » 1 . N 


66I2S 

IV* 



IF < I P 1 V < I >•» )Sw«76.S(, 


63 1 30 

2fc* 


SO 

00 7S J«1 »N 


66 133 

21 . 



IF ( IP | VI J»-l 155.7S.2Ss. 


601 36 

22 . 


ss 

1 NO * ( J » | 1 • N ♦ 1 


60 1 37 

23* 



IF t AWAx-OABSlM IND1 It <*v,»7S»7S 


GG1R2 

2 R • 


60 

I R" I 


C31R3 

2S. 



IC-J 


001 RR 

26* 



AMAX.OABSl A< l NO > > 


coirs 

27. 


75 

CONTINUE 
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I P I V C I C ) » I P I V C 1CIH 
I F M R*> I C >90,1)5,90 
90 S1GN.-S1GN 


33* IND*(L«I )*N*IR 

3 m * N»r c 

35 • AMAX»A(IND) 

36* r— A ( 1N0).A ( 1ND 2) — — 

37* 1 10 AUND2)*AMAX 

3 B* 115 INDEX (K,J ) *1R 

39* INDEX(K>2)»IC 

MC* r ND " UCM) »N*1C 

Ml* AMAX-AUND) 

' M2* DETErt-DE TER*ATf AX 

M3* JF (DETER) lMo.255, iMo 

~ M M * r M o AM NDl- l . 0D0 

M5* DO 150 U«l ,NN 

“ M6* 1 ND* l L- 1 ) »N«T C 

M7* 150 A( IND)*A(IND)/AMAX 

M 8 * Dtr 18 I t* I , N 

M9* IF (L-IC) 165, 181 , 165 

50* 165 IND*UC“-1)*N*L 

51* AMAX*A(IND) 

52* AT I ND ) *0*000 

53* DO 160 I*1»NN 

5M# I NTT* t I-T) *N*T. 

55* IND2- ( I-l ) *N* JC 

— 56* ATTN DT*A(ITJD y- ATTN D 2 T * AH AX- 

57* 160 CONTINUE 

56* 161 - CON T INU E 

59* 182 CONTINUE 

60* DO 235 1 * t, N 

61* U»N*1-I 

62* IR-INOEX tL ,1) 

63* I C ■ I NDE X ( L , 2 ) 


73* 
7 M • 
75* 


210 DO 230 K*1 ,N 

[ N^O*T T R ~T > • — 

IN02-UC-1 )*N-K 

AHAX-A(TND) 

All ND } * A ( l ND 2 ) 

230~ * TT N PZ) * A M AX 

,235 CONTINUE 

— D E TER ■ S IG N • DE T E f 

RETURN 

250" N**! — 

255 RETURN 


ENO-OT-T&MPI-tATK 


— NO- — DT*G NOSTtCS, 
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SJfORilS DhATh Ll UMMML 

MVJ tC 1 ?-. 9/ i6-|3:32 I #v | 


SUbHOUl inE OttATML ENTRY POINT 


STORAGE UStO: COOE,Jl CtS*?*? 0 *TA|u) OtuO S 3 ; bL*NK cOKMOl M 2) OLL-oCu 


external references (Block* name) 

Ga.03 IKaC*F 
tiCoM NERR3* 


STORAGE ASSUNNENt (BLOCK, TYPE, RELATIVE LOCATION, NaNE J 


EC J 1 

‘Jv v*w3 * 

)L 

Erci 

CLOl** 

i jog 

OOyl 

0 L 0 L 6 T 

I3SG 

ocyi 

toe in? 

I**3g 

0001 

L’OUl 35 

2L 

i.101 


3L 

.Lui 


9L 

L iOoC‘ 

0 OolCOO 

CO 

CiOD 1 

1 lCOgOS 

I A) 

OOtiO J 

LCOwCA 

IA2 

too- 1 

: :-l*wc2 

1 A3 

CCv’O I 

o::u-7 

I b I 

u-DlO 

1 C-uCClG 

1 b 2 

COCO 1 

1 COCO 13 

IB3 

OOlO 1 

1 W-0CLC3 

in 

Cjj« ] 

l 

IN 

wC uU 

C:^22 

HOP* 

J0i?0 

I CcOOlS 

LA 

o:c,o i 

[ uC'Oul* 

LB 

0000 ] 

i .DO. 12 

LC 

C‘30v 1 

1 U-CJll 

LM 

Cv*J« ] 

uce;i7 

LM 

„oco 

I OwOCI 4 * 

LP 








J 3 1 I 

U 


SUBROUTINE Dr, a TML ( C * A , b , M , N , k J 

jilt-1 

2* 

C* 

AbSl RACT 

v-K * 

3* 

C 

GENERAL MATRlA MULTIPLICATION ROUTINE V. J T H Th*NSPOS£ OPTIONS 

g-ici 

*»• 

c 

6HERE, M IS THE NUMBER OF ROjiS OF (a) 

v :i 3i 

S* 

c 

n is the number uf hoas of (b> 

Cuic i 

6* 

c 

r is ihe number of columns of t b > oR ibit 

; i i C » 

7* 

c 

TRANSPOSE OPTIONS A«£ CONTROLLED &Y ThE S I 6NS OF M AND N. 

vvIOI 

fa • 

c 

THE FOlLOMNG PRODUCTS may be obtained 

Oil 3* 

V* 

c 

IC)«(aMB) M AND N POSITIVE 

c-cici 

u» 

c 

{ c > - 4 A > 7 « a 7 M NEGATIVE FOR ( A ) T 

Jfll 0) 

i i • 

c 

(0«<AMB)T n NEGATIVE FOR ( fa ) T 

(.0101 

12* 

c 

(c»*«a»tib)t m and n negative 

1,0 1 u 1 

13* 

c 

*here t indicates transpose 

„ilti 

JM* 

c 

IF M is NEGATIVE, h IS The NUMBER of Rons OF (AJT 

CilOl 

lb* 

c 

‘IF N I S NEGATIVE, M IS THE NUMBER OF R 0 A S OF (0)T 

tfCICI 

16* 

c 


4.JI01 

17* 

c* 

OUTPUT ARGUMENT • C 

(i3IG3 

1 fa* 


C I MENSI ON C ( 1 1 

viioi 

19* 

c 


^^iu3 

2l • 

c» 

input arguments • a*b,m,n,k 


2 l • 


DIMENSION A t 1 ) ,B 1 1 I 

5-IC-H 

22* 

c 


CClO-4 

23« 

c • • * 


joins 

2M. 


DOUBLE PRECISION CD»C,A.b 

i'O 1 06 

2S* 


I *3»1 


?6. 


IH«IABS(MI 

;-»c 

2 7* 


JN*i ABSICI 

JCli I 

2B* 


IF(M , L T • C 1 GO TO | 

4,0 * i 3 

2V* 


I Al-IM 
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&04T*4 3L* 1-A2*4 


sans 

03116 

00 1 I 7 

31* 

32* 

lit 

1 

GO TO 2 

I A 1 ■ 1 - 

I A IN 

00 1 2D 
Cal 22 

3*4* 

35* 

2 

I F < N “,LT, o > GO TO 3 
I B 1 ■ | 

OOT2J 

36* 


fB 2 - I N 

091 2*4 

37* 


GO TO *4 

C012& 

38* 

3 

I B 1 -K ' 

00126 

39* 


I B2» 1 

03127 

*40* 

*4 

DO 7 LM<* 1 , I M 

00132 

*4 1 • 


LC-LM 

00133 

*42* 

. 

IB 3-f ' — 

001 3*4 

*43* 


DO 6 LP B l » K 

00137 

*4*4* 


CD ■ O.ODO - 

031 *»0 

*45* 


LA-I A3 

001M1 

*4 6* 


L B « I B 3 

001*42 

*47* 


DO 5 LN= 1 . IN 

aarir 

*4 8 • 


CD - CD '♦ A « LA ) *B { LB ) 

031*46 

<49* 


LA-LA + I A I 

001*47 

50* 

S 

lb-lb+ibi 

00151 

Si* 


C (LC J a CD 

00152 

52* 


LC • LC + IM 

00153 

S3* 

6 

IB3-IB3+IB2 

oats 5 

5*4* 

7 

— I A 3 * t A '3 * 1 A 2 

00157 

55* 


RETURN 

09160 

56* 


END - 


END OF COMPILATION; NO DIAGNOSTICS, 
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UMAP.IS EXPO 

Map »7M1-u9/1<,-13:33 


1. L I B SYSS«iiSFCS. 

2. I K* LOCEXP 


HsMqNI Tor ENTkY POINT TRACL 


ALkEADY 


ADDRESS LIMITS uOlOOO 037532 OljGuD 076 22o 

STARTING ADDRESS 031714 


WORDS DECIMAL 157 u 7 IBANK 


a ss i i dbank 


SEGMfcHT MAIN iiCUiOQ U3753 2 


NSWTCS/FOR 

I Du 1 Uliil 

0 y 1 D 2 I 


nrblks/msfcss 

l UU 1 1)22 

1 * 1 0 

3 

NRWNDS/F0R50 

1 DOl 1 1 1 

uu 1 1 7o 

2 

NWEFS/MSFC55 

1 Do 117 1 

Ou 1 4 2 B 

2 

NEXP6S/MSFCS7 

1 u j 1426 

001620 

2 

AL0GS/F0R51 

1 DO i 62 1 

0 J t 7 3 6 

2 

CSIG6V/SC402U 

1 uCi 1 7 3 7 

0G2225 

0 

*> 

CERMRK/SC4020 

1 002226 

002255 

4 

t 

y 

CLABLV/SC4u2u 

1 002256 

L- u- 3 3 l S 

£ 

u 

2 

2 

NFTCHS/F0R57 

t 003316 

0036 1 5 

NFTVS/FOR 

1 O 0 3 6 1 6 

Du 3 640 


NCL0SS/MSFC57 

1 003641 

004 UO 7 

2 

NWBLKS/MSFC57 

1 004010 

004 1 7 7 

■J 

N8SBLS/F0R 

1 004 200 

D04235 


nupdas/for 

nbfqcis/for 

1 Ou 4 2 3 6 

D 0 4 2 7 l 

2 

CTM0DV/SC4 o2D 

1 0042/2 

004317 

0 

-> 

CACCBY/SC4V/20 

1 u0432u 

Dd434 l 

4 

0 

•y 

CXM0DV/SC4u2o 

1 004342 

Ov4367 

4 

0 

■) 

CONCAT/MSFC 

1 004370 

d«454 1 

4 

0 

SET INT/SC4C20 

I 004542 

004574 

0 

5 

CH0LLV/SC4020 

1 uu4575 

004657 

4 

0 

2 


DEE I NED 


0*10000 076226 


G40U00 040301 
0*10002 010013 
010011 040034 
010035 0 4 0 1 0 6 
010107 040H7 
040150 D4C206 
BLANKSCUMMON 
010207 010222 
BL ankscommon 
0 4 0 2 2 3 C 4 0 3 3 1 
BLANKSCOMMON 
040332 040367 

040370 040H5 

040416 040420 


010421 042622 
042623 042631 
BLANKSCOMMOn 
042632 042642 
BLANKSCOMMON 
042643 042651 
BLANKSCOMMON 
042652 042673 
042674 042701 
BLANKSCOMMON 
042702 042716 
BLANKSCOMMON 


66 



CnONLN/SCm^Zu 


6046t>o 

00^2 6 5 

0 

;> 4 2 7 17 042770 





2 

bL ANKSCOMMON 

CL1NRV/SC4o20 

1 

ij|i ^266 

606127 

U 

042771 043070 


3 

GGG 


2 

BLANKSCOMMON 

rvQri u / «; r u r» ? »* 

i 

e l i U O 1— lit 

' - A 33? 

.. 

a.ll l tJ 7-1 tlt-1-1-1-7 

W-l-JVt-*-'— J W-I U * W 

3 

— uu u i — * u 

G G 6 

» J 4 

•J 

2 

J 1 •* w FI W f J * • t 

BLANKSCOHmON 

CXSCLV/SC462C 

1 

00*333 

306S35 

0 

04312c 043146 


3 

ggg 


2 

BLANKSCOMMON 

CERNL V/SC4020 

1 

06*536 

0067 1 2 

0 

043147 043161 





2 

BLANKSCOMMON 

CERRLN/SC402U 

J 

00*7 t 3 

00.7022 

0 

0-4 3 162 0431 77 





2 

BLANKSCOMMON 

CSETC V/SC4620 

1 

CO 7 0 2 3 

w 0 7 0 6 3 

rs 

V 

ij 4 3 2 0 0 04 32 1 j 





2 

BLANKSCOMMON 

N£XPSS/F0R57 

1 

00706*1 

007147 

2 

043211 043220 

CSETMV/SC4G2G 

1 

007 1 5jJ 

007226 

0 

043221 043235 





2 

BLANKSCOMMON 

CFRAM/SC4620 

t 

007227 

007446 

3 

043236 043322 


3 

GGG 


2 

BLANKSCOMMON 

C X A X 1 S/SC 402u 

1 

007447 

0 j 7 *52 

0 

043323 043356 





2 

BLANKSCOMMON 

VCHAF'V/SC4y2U 

1 

O07 6S3 

.10*16 

0 

0-43357 o43376 

RlTE2v/SCHj2u 

1 

010117 

0 10347 

6 

•J43377 643425 

BPL0.TK/SC4u2L 

1 

0 1 o 3 S 0 

010^30 

6 

043426 044065 





2 

BLANKSCOMMON 

CCAMRa/SC4u2G 

i 

0 1 u 4 3 l 

0lo&02 

0 

044006 044G17 


3 

GGG 


2 

BLANKSCOMMON 

TABlI V/SCHw20 




0 

044020 044340 

NBDCVS/F OR 57 

1 

0 1 0 5o 3 

010*36 

2 

U 4 4 3 4 1 0444(10 

N C N V T S / F 0 R 5 7 

1 

010637 

j 1 1 o 7 | 

2 

j 4 4 4 0 1 344470 

NOT I NS/MSFC5B 

J 

ui icy 2 

01 1422 

2 

■34447 1 644501 

N0UTS/F0R57 

1 

01 1423 

<. 12377 

2 

0-44 50 2 044532 

NI0ERS/MSFC57 

1 

01 2400 

012557 

2 

644533 044655 

NININS/MSFCSS 

1 

U 1 2560 

C 1 3ol 1 

2 

u 4 4 65 6 64467 7 

NI NPTS/F 0RS7 

1 

013012 

013671 

2 

044700 044722 

NFMTS/F0RS7 

1 

013672 

014576 

2 

044723 044741 

NFCHKS/MSFCS7 

I 

014577 

015415 

2 

644742 645116 





4 

045117 045176 

NT ABS/MSF CS5 




2 

045171 045257 

CPLOT V/SC402w 

l 

01 54 i 6 

0 15567 

0 

04526U 045317 


3 

ggg 


2 

BLANKSCOMMON 

CLINEv/SCHi/2u 

1 

015570 

Jl 6 1 S3 

0 

■045320 045376 





2 

BLANKSCOMMON 

TSCLV 1 /SC4t,2G 

1 

Ol*i54 

i*l 6272 

0 

645377 0454Q7 


3 

GGG 


2 

BLANKSCOMMON 

XSCLV 1 /SC4G2G 

1 

01*273 

C 1 6 4 1 r 

0 

C 4 5 4 1 c 045420 


3 

GGG 


2 

BLANKSCOMMON 

CAPLCT / 5 C 4 o 2 b 

l 

01*411 

61**53 

Ci 

04542) 045474 


3 

GGG 


2 

BLANKSCOMMON 

CAPRNV/SC4620 

1 

01*654 

01*736 

0 

045475 645513 





2 

BLANKSCOMMON 

CPRNT V/SC4L2G 

1 

016737 

Ol733o 

0 

u 4 5 5 1 4 0 4 5 5 5 y 


3 

ggg 


2 

BLANKSCOMMON 

CGKDl V/SC4G20 

1 

017331 

020235 

6 

0 4 555 1 645*64 


3 

GGG 


2 

BLANKSCOMMON 

CDXDYV/SC4o2b 

1 

020236 

02 1 03 Q 

0 

b 4 5 6 6 5 045/51 
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SYS$*RLIB$» LEVEL M57-0 

ENO OF COLLECTION - Time 3.972 SECONDS 
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-SXQT -E.XE-Q- 
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Typical Output Results for Set 1 
Data Using a Two-Term Exponential Model 



• OOw C* GO 01-006 C'i;0GG. - G0 + 0iOU» 

» ^CcOji’ jCLiuJDGui'.'jiii-O + Qjl t 

.79999999999999V999D* C G) , 

• 12 c i C 0 J S> (iiroZZ (j l. j iTC 0 >' QC‘ 27 

• lfo»wjOCjUtjUwuC'v’^C'G^ + CL*2 1 

• 2 0 a 0 uct I- a j c c e c :• -j o c- c * g 0 2 , 
. 2«iaes-3:c«aoocGix-c B *oc2, 

• 2 8 :j 5 0 J £ C o G w w 0 C 3 G O + 5 j 2 , 

• 32*}COO0O0OoOvOCGvGO + O32 , 

• 3 6 w t GoOOOO 6 Gi>0 « 60 ♦ C u 2 f 

• ‘iOcOaCCOGGuOCGui;w : iiO*oa2 , 
. 9399999999999 999990* u< j2, 

• 485000 0C.OCCO 6- g6G*->00 + C02 , 

• 5 2 Vw(-Cuuui uwGv'Lilj + L;u2 • 

• 5 6 w t *j j L- T u C w' 0 C t ii J >.j w P * t. 1 c 2 t 
.5999999999999999990*002, 
.63 9999999999 9 9 99990* GC 2, 
.6799999999999999990*002. 

• 71999999999999999 9 D + uii 2, 

. 7599999999999999990 * 2 1 
. 7999999999999999990* j&2 , 
. 83999999999999 99990+JJ2. 
.8799999999999999990+052, 
,9|99999999999999990*oo2, 
,9599999999999999990*052, 

. 10'jCCGGwGGwJwGGG B i50 + 09 3, 

, 1 0*0C00CG0j0C'GGOGC0*u03 , 

• 1 3 8 5 y O u 0 0 0 U u i. 6 c uO 0 + G ■- 3 , 
• I l2cc0OC’C0UvOGOGCU0 + CiO'3, 
« 1 t66C0GC-uG*+GwCwGO0D + C03, 
. I 2vCG^vCCC<ilbbvuvvG^ + Gu3 , 
« 12*l50v-^'00Lo05ii5'c»’j0 + 5-j3| 

• 1 28U57^!'UCuiItiOv2vu^*L>3 , 

• 1 32goGCUvcwUC 6 *, c i j 0 + 6 l’ 3 , 
. CGOCOOOOOOCwOCvOGOP+OOO , 

• OCX OCOOOoCoCfGCivJ- 00+006 * 
. v CC'OQCCC GOOO i Go 6wG0 + G v 0 , 

• CC.X'CO+GOOuOcCCG-^wO + C'Oo, 
. SCOGDCiOOOOwCOCOCOC'O + c-SC, 
. o 0 G L< 0 ; 3 O 6 GuOucGC 6600+0 06 , 

.0&SCCC8t*»C«OtOOe5CO*CCO, 

. 0 0 0 (• 6 J V 6 0 1 6 G‘ G o 6 u , 0 0 + O 0 u , 

• GO JCO€XGOGGCXC?C'i7-J00 + O3G , 
. o 0 _* 0 0*0 OOu G ,J,Ob Ou uO + oOG , 
. GCGOCOOOGGGOCOGOOC-P + GOG , 
. wDoiiOO jCwGuOtGO.’v CO + OCC- • 

• CPOCGC’OCuOOt'OC-COOOO + CJu , 

« »«0'i0fcCG0uG60oGv J ,tD + 6 00 , 

• eOOOCOOC (JCOOOt; 00 jC'O + G 00 , 
. 11 0 C 0 0 b e 6 0 fi 0 W L J i; 6- u 0 ♦ 0 J 0 , 
. 1 3 jtcacaooGaoOoL-^i £>*oo i . 

.61 OCOOjOiiuuOuOjCOt*^ + COO, 
.9 70o0a'ul?U£J000CC000 0 *C30 . 
. 305(>C'!i3t J0O5>i,c0 j^u^ + OOO , 


,2c>OCoooCuOoCl>OOOOOD + OOI i 
.5999999999999999990*001 1 
. I OOOOOOCCOO SCOv C OOP + 00 2 1 
.1900005 50 OoC 0 C OC 00 0 ♦ 002 , 
. 1 BOOOCGCGCiCOCGCOOO+OCZ , 
.22oeoutoocoeocccooo*oo2 1 
.260000CcCCuCC'Co0°CD + 002 i 
. 3000oCC0CC0GC0uClCCD + 002 1 
.3 400000000000000000 *002 
. 3eOOG0CCOC„C0tuOOOD+OO2 , 
.9 200COOOUCOOOODCCOD + OC12 , 
. 4600006000600 OOC 00 0+002. 
,5000000 C 00 000000000+002 
.5399999999999999990*002 
,579999999999999999d+002i 
.6199999999999999990*002 
, 6 60000000000 JuoCOOO* 002 
.7UOOOOCOCCOOOOCCOOO+O02i 

. 7400000COCGOOOCOOOD+002 , 
. 7BOOOOOOCCillioOuCOOO*002 
.82oooaocootioococooD*oo2 
.8600000000000000000*002 
.90000C0C0C0C0C0C00D+002 
• 9400000000600006060+002 
. 98O00000000C00CC0 oD* 002 
. 1 020000000000000000+003 
. I 0&0O0003CO000U0O00+0C3 
. * lOOOOOOOOOOUOJOOOD+003 
. 1 I4COOOCOOOCOOOOOOO+003 
. 1 1 B000o000uC0t0000D*003 
. 1 22COCOOOOOOCOOOOOD+003 
. 1 Z 6 OOOOOOOOC OUOOOOD* 00 3 
. 1 3OO0CC00C0CO0CDOOD+003 
. I 3 4000000 00 COO 00 000 + 00 3 
.DOODOOOOOCCOPCOCOGD+OOO 
.OOOOOOOOCOOCOOOOOOD+OOO 
.GGGCuGGOOOoCuO 00000*000 
.COOOOOOOOOuCCOUOOOD+COO 
.OUOOOGOOOOGCOOOOOCD+ODG 
.0O0C0O00000C 0000000*000 
.OUOOOOOOOCOCOCCCOOD+OOO 
• OOOOOOC OOOO'OOOCCOOD + 000 

. ooooooooooocooocooo+ooo 

, GOOUOOUOGGOC 0000000* 000 
, 000000 OOOOuOCCO&OOD* 00c 
. OOOC0OOOOOOOCO0OOOD+OOO 
.0000000000000000000*000 
. OOOOOOOOOOuO 00 00000*000 
.COOOOOOOCCOCCOOCOOD+COO 
. COOCOO0OOO0C OOOCCOO + OCiO 

. 7400o0c(j00oc-00000c0*000 
. 53oOOOOOOCO0OO>aOOoD*oU O 

.42SOOOOO0OOCCCOGOCD+OGG 
. 345COOOG5C'OrOCXoCOD*OoO 



3 i 

25^00 CwOOi>Ct'uO^:>wC»0 + GOo 
2 ! ^OoOwCtUuvuiiC-w jU^ 4 ll>u 

I 85C'wCUNi)w0uoi#i.u»j^ * t sj 6 
1650v0iCCCJCi.Co^«CD + 0CC 

i 0 w‘ (.• C ki w L> c w c? c* v w D ♦ £ 0 i> 

i H UCCCC&uOwutf Ofc-JCU+CCO 

) vu j^CuJ 

1 25Q£»£C»Q^Q6Ji4/Gfcii?*«jo^ + C36 

i i 5CCOC3oC^^wCytyv D+ oJG 
i 1 -jOUvil6ot<Joi,cvvV^ + iJl/U 

i o^u^rccucoc OGCCwCP+CGu 

1 0 wC ot/ji'Uui, Ct/GOG w 0^ + Oiiu 
1 Q .. CO t! 0 00 C wCft 30C J» v ^ *C C o 
9499999999999999990— ftft 1 

9M9999999999999999O-00I 
8999999999999999990-CuI 
8999999999999999990 -Cu I 
8999999999999999990-JOi 
8999999999999999990 -OjJ 
7999999999999999990-^^j 
7999999999999999990 -ft „J 
7999999999999999990 -„J| 
799999999999 9999990-ftft | 
7999999999999999990-ft ft | 
79999999999y999999D-ftft| 

7 4999999999999 99990- uU l 
749999999999999V99D-GG1 
7499999999999999990- w ol , 
iiOcaCGOOWCCCvCGCO'.O + JCC 

wOOGGOOvOCOOuGCCCOP+COC 
^ 0 OC G 0 w ^ G 0 — 0 ^ u * j ■* -j 0 ♦ o 0 6 
uO'j 00 COC GO CCCCOw GO 0 + GOO i 
0 0 m G £ G — 0 tl 0 J C v 0 o 0 o o 0 + 0 G 0 i 
■iD^C JO^icUuw • GO o C wCo'0 + ftOj i 
wG-i)C{>0oOoj0uOO0ut.0'fft , Jft j 
oC-.'v£j jO-Jvio oOwCOO^.C-0 + C.C0 i 
uOi.OC^OOCOsiO*CA>CCuoO + ftOOi 
C0oG00G3G0vCc 0 jo v * Jo j i 
uCmGoGO'jOOwJoC 1> o o 1/ 0 + ftuC i 
-COCwGOCOuOliOvGOCOC^*uOU j 
wG- GPC oCCOwu t C v OC i 

1 w370y^uvuwut c G w 0 0 ♦ ft v ft i 
» ubf^bCh yiiyOLjfOwl 0^^ + UoO i 


NN 

c 

♦ 68 

TL 

= 

• ^ o C i# u £ u 0 E ♦ V* 0 

TR 

= 

• lbCG'C3G0E + ;3 

Y 8 

c 

• « w 0 U w 0 it 0 E ♦ G i‘ 

YT 

= 

• 2 t 000 jOrE + G | 

Yb 1 

= 

- . 1 c 0 C w . u 0 E ♦ ft j 

YT 1 

c 

• icwv>w: : C0£ 4 C 1 

VARY 

s 

• ^ wL Jw »-v0^ ^ i 1 

TOLER 

C 

• li.wo ! u s uuE + C5 

I T E R H 4 

= 

♦0 


.2aoCoocaoOGOt)CGOOCD + OOG 
. 2 JGOGOGGoOCGOGGOOOD + eaC 
. I 950CCC!QUG<j0CCJCl!Q0D + 0 30 
. I TBOuOCGOGoC CocGOOD+OUO 
.IbSGGOGOC’C'GOGCGOOCD + GGG 
. 1 4bCCOOOOGuOuOaGOOo + GGO 
. • 3bCGJOOOCG£IGCCOauo + OGG 
, I 2S00GC0vJ0o00C/OLiU0D + 0 JO 

• *20GG0C000oC00uo0uD + 000 
. I 1 BOftoftGOOuGuOOGOUO + CGO 
. I loGOOOOOOjCGGoOaoO+OGG 

• l05GwyU0uOut0o0w0u0+tijjy 

. * O&0OUG00G0CC0ft0Ca0 + 000 

. ICOOCOOCoOoCoCCJOOO+CftC 

. *U0GftJGCo0oC0DGC0u0+0C0 

• 9899999999999999990-GO I 
.8999999999999999990-001 
.8999999999999999990-001 
.8499999999999999990-001 
.8499999999999999990-001 , 
.8499999999999999990-001 
.7999999999999999990-001 
.7999999999999999990-001 
.7999999999999999990-001 
.7999999999999999990-001 
.7999999999999999990-001 
.7499999999999999990-001 
.7499999999999999990-001 
.7499999999999999990-001 
.7499999999999999990-001 
.CU 000000 OGOCOOOOOOO+OOO 
.OUOOOOoOOOOtii' 0^000 0 + 000 
.OOOOGOCOOOul-OOOCOOO + ooC 
.OOUOuOoOOOOGOCCJCiOD + OOO 
. CCOCOGoGOOGCjCOCOOD+GCO 
.COilOoOCOOO^COGOOCOO + OCO 
.cuooGaocGoocGcoooaD+ooo 
•GOUCoOGOCOjCUCcCOGD+GjU 
« GCOCuCOC OCCOuUbCCcD + 000 
.OOoCuCoOCCoaOUGOOuO+OOO 

.OOUCoOoGOGuOOGuG'OOO + ggo 

• OGOCOOCC OcjCCCoOOGD + goo 

.GOOPOOGCCGGGOCJGoOD + OOCI 
.OOOOOCoOOGwCGGGGCG 0+000 

. CG10C0 0 00 00 01- GO OG 00 0 + ft 0 0 

. GOOD OOOOOv oC'GGc COO 0 + COO 
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INITIAL ESTlMATts F0« T'60 -tEkM EXPONENTIAL MODEL 


I 

F 

A 1 

8 1 

A 2 

B 2 

I 

.590963*^1 

.215863+01 

. &C l 193-01 

. 75o0C0*6 1 

-. 1 98030- 1 6 

2 

.19671 9*U 1 

• 152969 + ol 

. 7 1 6359-0 * 

.750600-01 

•ocoooo 

3 

. 5C-C553*i.L 

. 12,9223 + 01 

.626321-01 

. 7 5 u 6 C 3 — 6* 1 

• ccococ 

9 

. 339888 + ,,,. 

.6877 39 + 0-0 

.529697-0* 

. 7500CG-C1 

. c C OOoC 

5 

• 7 5 1 ; 3 S ♦ Q 

. 9368 | 5+uo 

• 910629-C.l 

•756006-01 

-.362522-1 6 

6 

. 1 779 15 + i,l 

. | 9 39 9 7 + 0 j 

•262821-C* 

. 1 66386+60 

.26891 1-02 

7 

•191917+ti 

. 1 CO 989 + 0=3 

. 1 652 19-Cl 

. I 22265 + CC 

.376976-02 

B 

. 1 9L5u9 + 0 1 

•978 v 26- t 1 

.197689-01 

. 127769+60 

.910700-02 

9 

. 1 8 6 9 2 9 ♦ i- 1 

• 1 2 1 1 7 8 + 00 

. 1 56778-0* 

• 1 2 7 7 6 9 + 60 

. 9 10700-02 

1C 

. i92<,7 1 +01 

•971 186-01 

. 1 73739-01 

. 125225 + 0 0 

.399909-02 

1 1 

. 1 8 1897+o 1 

• | 18973+00 

.211909-01 

. 121696 + OU 

.372338-02 

1 2 

. 179397+01 

•136673+00 

.25 w 953-0* 

. 1 1 7999+60 

.397515-02 

1 3 

. 166707+01 

• 1 563 39 + 00 

.299937-01 

. 1 1 9 3 3 6 + 60 

.322693-02 

1 9 

. 152889+01 

.1892 I 9 + „o 

•3922C6-0* 

. 1 1 1CC3+C0 

. 29B965-02 

1 5 

. 1 37230 + 01 

.2285 18+00 

.390307-0* 

. 1L7997+G0 

.276820-02 

16 

. i 2 2 o 9 9 ♦ u 1 

•272315+Ou 

.959760-0* 

. 10531 9 + l-C 

.256922-02 

1 7 

• 128961 + 01 

. 298992+Lo 

.929^99-01 

. 1 C- 9 1 56 + Cil 

.285635-02 

I 6 

. 1292u 9+ w | 

. 299959+ ue 

• ‘*20069-0 * 

.111971+06 

.302638-02 

IV 

. 126053 + 0 1 

. 25 1 692+OC 

.929768-0* 

. 1 1 2668+00 

.3 1 1690-02 

20 

. 1202 o 5+01 

• 266993+uo 

.951299-0* 

. 1 1 3065+60 

.319559-02 

21 

. 1 1 1593+0 1 

• 29 1 365 + 00 

.986171-0* 

. 1 12902+00 

.313353-02 

22 

. 1 3L So8 + C 1 

• 329919 + oO 

.515887-31 

. 1 12356 + 60 

.309299-02 

23 

.105693+01 

. 3 1 37 9 3 + Uo 

.591731-0* 

.115 1 05 + 00 

.329661-02 

29 

. 1 utb63+o 1 

•3199o2+00 

. 527675-0* 

.11691 l+oo 

.392868-02 

2 5 

.9&8l39 + i.f 

•3351 19+oC 

.598990-0* 

. 1 1 7983+00 

.350610-02 

26 

• 989789 + c-L 

.322318 + 0-0 

.596993-01 

. 121272+00 

.379201-02 

27 

. 98 C 938+00 

•322770+00 

.555779-01 

. 123609 + 06 

.390656-02 

28 

. 969877 + OL 

• 328668 + oC 

. 597796-01 

. 1 25 166 + 60 

.90151 1-02 

29 

. 930089+00 

.333333+00 

.569731-0* 

. 128359 + 60 

.923981 -02 

30 

.9^8669+00 

• 3 38991+0-0 

.589712-0* 

. 1 36677+60 

.939175-02 

31 

.899165+00 

• 3569 12 + 00 

. 6028 1 9-0 * 

.1 32265+oc 

.999839-02 

32 

•779io9+oO 

•382289+ou 

.637355-0* 

. 1 33255 + OC- 

.956952-62 

33 

. 75L235 + 0C 

.392667+00 

.666289-0* 

. 1 365 1 1 +C0 

.971982-62 

39 

. 683900+«L- 

•9171 C 3+ vO 

.689371-0* 

. 1 37 105 + L6 

.961995-02 

35 

. 6o5232+ w o 

.951929+60 

.732927-0* 

• 1 3 8 1 9 5 ♦ 6 0 

.988777-02 

36 

. 55552 0+wo 

•977088+00 

,. 792927-01 

. 190173 + 60 

.502029-62 

37 

. 976859 + uc 

.51 2865+06 

.8197 1 6-0* 

. 19 16 18 + 00 

.51 190&-C2 

38 

. 92&966+u0 

.S9‘.-759 + t-6 

.857185-0* 

. 1 93836 + Lk 

.525707-02 

39 

. 399652+00 

• 59 237 9 + i-y 

.927223-d 

. | 95973 + 00 

.536189-02 

9 Cl 

. 3C996C+0L 

•62189 l+o 0 

.952586-0* 

. 1 97729 + 00 

.550622-02 

91 

. 22® 0' 19 + Oo 

. 69C259+0D 

• 1 039 80 + 00 

. 1 99922+60 

.661205-02 

92 

. 1 8“ 256+oL 

.739895+10 

• 1G7960+CU 

. 161626+00 

.579981 -02 

93 

. 1 33032+01 

• 829925 + 60- 

. i 1538 1+Ot 

.153262+06 

.565332-02 

99 

.953392-01 

• 921 1 0 9 ♦ 0 o' 

.127393+06 

. 1 56362*01 

.698209-02 

95 

.92856C-OI 

. 9299S2+vC 

• 126726 + QC 

. 157760+00 

.6 1 2939-G2 

96 

.9o5l9l-„| 

•975897+00 

. * 31825 + Cj 

•160397+06 

.628977-02 

97 

• 10689o+0L 

• 1 3 7 C 5 0 ♦ 6 1 

. 1 92296+CL 

. 163706+60 

.695887-02 

98 

. 15991 6+oL 

• 1 182 19 + v 1 

. 153529+06 

.1661 33+CC 

.663315-02 

99 

. 19L773+CC 

• 1 15266 + 61 

. 152772+01 

. 169139+00 

.680981 -02 

50 

. 156lo9+oo 

. 1 18221+01 

. 158558 + Q6 

. 1 72752 + 60 

.702008-02 
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Si- 

52 

—.697.7.93^64— 

.827253-01 

— +1-62191+01 
.951 721*00 

• 1.45658 +00 

• t 4Q500 + 00 

--.-1768 47*00- 
.181307+00 

.725439-02 

.750509-02 

S3 

.681064*01 

.886495+00 

•134412+00 

. 1 86036+00 

.776607-02 

54 

.980894-01 

.840089+00 

. I 30080 + 00 

.191381 +C0 

. 8055 13-02 

55 

C A 

.106007*00 

_ I nAliCt 1 4 A n 

.810961+00 

. T Q -» Q T O a ■> n . 

. >28233+00 

. 197582+00 

.838273-02 

3 O 

— • 1 ^ofa-i-’OC — 

. 7 4 7 87 4 ♦ U0 

• • 2+563+00 

. 209789+00 

.875336-02 

57 

08707*00 - 

+790074+00 

“ • 1 3233 6+00 

•212778+00 

.915148-02 

58 

• lOS861*oo 

. 789423+0C 

. 1 37549+00 

. 22 1 65o*oO 

.958047-02 

59 

. 995807*01 

* 798074 + 00 

• 1 46680 + 00 

.231457+00 

. 100378-01 

60 

.906543-01 

.819589+00 

. >62410 + 00 

.242216+00 

• IC52I2-C1 

6 1 

.838456*01 

.850950+00 

. 1 85905*00 

.253707+00 

. 110178-01 

62 

.840578-01 

•89o’l8+00 

• 220787-fOO 

.265941+00 

• 1 15260-01 

63 

. 1 0061 T *0C~ 

. 944 &cl* 00 - 

~ • 2 8T'7 30*0tfr” ” 

• 2 7 9 24 5+00 

• 120566-01 


The initial estimates USED are 

■ S 2 
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NON-LINEaR EXPONENTIAL REGRESSION ANALYSIS USING an ITERATIVE 


ASSUMED model IS T«0 EXPONENTIALS 

A I * .951721 ♦(/(> 6 I c .I90SCC*00 

observed response, computed- response , 


A2« ,I6I3cT + 0C B 2 a 

and RESIDUALS Using INIl 


QB SfcR V£p_ 

1 30000*0 1 
790000*00 


rnMOIITcn 

v V “ ' W T CTV 


1 l33(,3*Ct 


610000*00 
530000 *00 
970000*00 
925000*00 


6971 82*00 
718990*00 
582961*00 
980029*00 


RE-5-1 DUAL — 
. 166972*^,0 
•i 15 7 | 8 2 ♦ 00 
-. 108990*00 
-.529608-0! 
100290-01 


> 385000*00 

• 395UUO+00 
.310000*00 
•280000*00 
. 2 55000* vO 
•230000*00 
.210000*00 

• 1 95000*00 ' 

• 185000*00 

• 175000*00 
•165000*08 


392007*00 
296396*00 
261303*00 
239285*00 
213335*00 
796973* 00 
189086*00 
173828*00 
165569*00 
158813*00 
153213*00 


.929933-01 
.98659C-C1 
.986970-01 
.657 I5C-01 
.916659-01 
.330255*01 
.259195-01 
.211721-01 
. 199359-01 
. 161866-01 
. 1 17873-01 


5000*00 — 
150000*00 
1 95000*00 
190000*00 
135000*00 
130000*00 

125000*00 

125000*00 

120000*00 


199932*00 
1 90888*00 
137738*00 
139892*00 
132283*00 


.556829-02 
.91 1223-02 
.226237-02 
.108072-03 
-.228305-02 


127589*00 -.2S83S7-02 

125925*00 -.592989-02 


120000*00 
1 I 5000*00 
1 15000*00 
m>eo6+Do^ 


. 1 10000*00 
. 105000*00 
. 105000+00 
. 105000*00 
. 100000*00 
• 100000*00 


• 100000 + 00 
. 100000*00 
. V50000-01 
.950000-01 
•950000-01 
r5CCOCO-& t - 


123362*00 
121377*00 
1 19957*00 
I t 7599*00 
1 15779*00 
11 9006*00 
1 12271*00 
I 10572 + 00 
108909*00 


-.336156-02 
-.637667-02 
-.995738-52 
-.759387-02 
-.577866-02 
-.900600-02 
-.727195-02 
-.557158-02 
-.890371-02 
-T726S 72-02 


105656*00 -.565599-02 
109073*00 - • 907 302-02 
102516*00 -.751589-02 
100989*00 “ -.598351-02 
999752-01 -.997529-02 


•900000-01 


965283-01 -.652833-02 


CORRECTION procedure 


.750509-02 
lAU ESTIMATES 
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• 9U0uQU=CU — 

«-9 SG884-C1 

508860-02 

• 850000-01 

.93870-7-01 

— .867072— q 2 

•850000-01 

.922793-01 

-.727927-0 2 

.850000-01 

. 908989-01 

-.589886-02 

• SSqCo G^U l 

.8959*4 1-0 1 

- .95991 1 - 02 

•850000-01 

.882097-01 

-. 320967-02 

•85eooe-6t 

r 8 68 9 5 ?— 0 1 — 

-» 1 B95 2C-0* 

.800000-01 

. 856009-01 

- • 560C9 1 -C2 

•eoeooo-oi 

. 893250-01 

- .932995-o2 

•BucCGC-oi 

• 8 30685-0 1 

-.306859-02 

.socooc-oi 

-.818309-01' 

— . | 83c88-q2 

.800000-01 

• 8y6 1 17-01 

— .61 1 679-03 

• i ' 

- .799-tar-o i — 

— *589337-03 

.800600-01 

.782276-01 

. 177295-02 

• 8«0ii00-ii 1 

.770*21-01 

.293792-02 

.80Cv00-Gl 

.759190-01 

.908602-02 

. 8UCU0O-U1 

.797830-01 “ 

.521701-02 

.800000-01 

.736688-01 

.6331 15-02 

. 780000-01 

.7257rr-0T~ 

•".292869-02 

• 750c00-Ul 

.719901-01 

.350987-02 

. 7 5y OCO-0 1 

.70925 1 -0 1 

.967999-02 

• 75C606-ul 

.693759-* 1 

.562919-02 

« 7 SQOCC-0 1 

.683923-Ct 

.665771 -02 

* 750000-01 

.673291-01 

.767589-02 

- . 75CO00-01 

- .66321 | -Cl 

.867889-02 


SIGYIN. .359S25-C1 


CYCLE NUMBER » I 
B matrix 


. lOCOGG+Ol 

. GoOOCO 

. 1 CGC06+C 1 

•000000 

• 755028 + 00 

-.1937 15*01 ~ 

• 985102 + 00 

-. 35721 1+CO 

.570068+00 

-.217018+01 

. 970926 + qo 

-.703779+00 

. 9309 17+yO 

-.295782+01 

,955968+00 

-. 1Q3999 + 01 

. 329977+00 

-.297930+01 

.991726+00 

-. 136593 + 01 

. 295367+uO 

— .23 352 1 +!> 1 

. 927696+00 

— * 168198 + 01 

. 185259+yQ 

-.211578+fcl 

. 9 1 3875+00 

-. 198830 + Cl 

• 1 3 9 8 7 6 + GO 

— « 1&6372+01 

. VC026C+0S 

-•228513+C1 

. 1*5610+06 

- . 1 608 l 8 + 0 1 

, 886898 + 00 

-•257266+01 

. 797386-01 

-. 1 36 60C+O l 

. 8736 36 + 00 

-•285113+01 

. 602o99-0 1 

-. 1 1 9597 + 01 

. 860620+00 

-•312073+01 

.959569-01 ‘ 

-.951761+00 

.897799+o 0 

««338| 66+0 1 

. 3932C?-ul 

-.783939+00 

.835168+00 

-.363912+01 

.259132-01 

-.691217+00 

.822726+00 

-.387831+01 

. 1 95652-01 

-.521378+00 

. 8 1 0969+00 

-.9 1 1 992+0 1 

. 1 97723-y 1 

-.921773+00 

,798399+00 

-.939263+01 

. 1 1 1535-01 

-.339681+00 

, 7 B 65GC+0C 

-•956313+tl 

.892121-02 

-.272998+00 

.779782+00 

-.977609+01 

.635825-02 

-.217896+00 

. 763290+00 

-.998170+01 

. 98^066 - c* 2 

-. 1736J8 + 0G 

.751869+00 

-.518012+01 

. 362963-0 2 

-. 137986+00 

.790667+00 

-•537152+01 

. 273670-02 

-. 109392+00 

.729633+00 

-.555607+01 

. 2*6629-02 

-.865273-01 

. 7 18763 + 00 

-.573393+01 

. 15601 0-o2 

-.683301-2! 

. 7p8o55 + oo 

-•59C525+C1 

. 1 17792-02 

-.538106-01 

.69750 6+oo 

- >607020 + 01 



t_8 8‘9 36S-Q3 

-.92321 9-C 1 

.687(1 9.QQ 

- .627892*01- 

.67|996-03 

-.332320-01 

.676878*00 

-.638157+01 

• 506998-03 

-.260561-01 - 

. 666799+00 

-•652828+01 

.382798-03 

-.209018-01 

.656860*00 

-*666921+01 

.289023-03 

-.159590-01 

.697079*00 

-.680999+01 

.218221-03 

-.12961 |-01 

.637933+00 

-.693926+01 

.^<H*763-tr3 

-.972212-02 

— r6r27-9-3T’f 00~ 

--.705865+01 

. 129901-03 

-. 757726-02 

. 6 18582+00 

-•7 17780+01 

.939260-09, 

-.559983-02- - 

.609366+00 

-.729182+01 

. 709 168-09 

-.958952-02 

• 600288+00 

-.790086+01 

.535992-09 

-.35671 9-02 

.591395+00 

-*750503+01 

• 909279-09 

-.277029-02 

.582535+00 

-.760996+01 

. 3CS-2J8-09 

-. 219? 7 1-02 

" i i573 8 56-+ oo — 

**769925+01 

.230963-09 

-. 166696-02 

.565307+00 

-.778959+01 

. 1 790C6-C9 

-• 129172-C2 

•556885+qq 

-•787592+01 

. 131380-09 

-• 100029-02 

.598588+00 

-.795702+01 

.991959-05 

-.779132-03 

.5909 15+00 

-.803999+01 

.798963-05 

-.598798-03 

.532369+00 

-.810778+01 

- • *56&9fr 1 — u5 

-.96 28 3S -03 

- -r5 79933*cnr~ 

-.817716+01 

. 926959-05 

-. 35758C-03 

. 5 16620+00 

-•829267+01 

•322362-05 

-•2761 19-03 

.508923+00 

-•830991+01 

.293393-u5 

-.2131 H-03 

.501391+00 

-.636298+01 

. 18 3 76 8 -05 

-. 1699C2-C3 

. 993872 + 00 

-.891698+01 

• 1 38750-05 

-. 126770-03 

. 986519+00 

-.896800+01 

• 109 7 60-05 

-.977286-09 

.979266+00 

- -.851563+01 . 

.790970-06 

- - • 752783-09 

.972 1 26 + 00 

-.855997+01 

♦597205-06 

— .5 797 90-09 

. 965092+00 

-•860109+01 

. 95Q9Q7-06 

-.996303-09 

,958163+00 

-.863908+01 

.390997-06 

-*393952-08 

.951338+00 

-.867909+01 

.257097-06 

-•269208-09 

,99961 9 + OC* 

-.870603+01 

» 1 99078-06 

— . 2o3» 79-09 

*937990+00 

-.873515+01 

. 196539-06 

-. 156195-09 

.931969+00 

-.876197+01 

.11 0638-u6 

-♦ 1 20038-09 

.925036+00 

-.878507+01 

.835395-07 

-.922218-05 

. 9 1 B7uR + uO 

- • 88Q6C1+01 

.630709-07 

- .708306-05 

, 912966+00 

-.882939+01 

.976203-07 

-.593855-05 

.906321+00 

-.889026+01 

. 359597-07 — 

-.917570^05 

.900268+00" 

-.885370+01 

.271968-07 

-.320369-05 

. 399305+00 

-.886977+01 

.209966-07 
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PARAMETER corrections 
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TOLER* -.126392-03 
END OF CYCLE NUMbER 1 
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